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fi 要 结构 方程 模型 是 心理 学 、 管 理学 、 社 会 学 等 学 科 中 重要 的 统计 工具 之 一 。 然 而 , 大 
时 使 用 结构 方程 模型 的 研究 在 研究 设计 和 结果 报告 的 环节 忽视 了 对 该 方法 的 统计 检验 力 
行 必要 的 分 析 和 报告 , 在 一 定 程度 上 降低 了 这 些 研究 的 结果 的 证 明 效 力 。 结 构 方 程 模型 的 统 
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计 检 验 力 分 析 方 法 主要 有 Satorra-Saris Y% ~ MacCallum 法 与 Monte Carlo 法 三 类 。 其 中 Satorra- 
Saris 法 适用 于 备 择 模型 清晰 、 检 验 对 象 相对 简单 、 检 验方 法 基于 分布 的 情形 ; MacCallum 
法 适用 于 基于 ?分 布 的 模型 拟 合 检验 且 备 择 模 型 不 明 的 情形 ; Monte Carlo 法 适用 于 检验 对 
象 相对 复杂 、 采 用 模拟 或 重 抽样 方法 进行 检验 的 情形 。 在 实际 应 用 中 , 研究 者 应 当 首先 判断 
检验 的 目的 、 方 法 以 及 是 否 有 明确 的 备 择 模 型 ， 并 根据 这 些 信息 选择 具体 的 分 析 方 法 。 
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1. 引言 

结构 方程 模型 (Structural Equation Model, SEM) 是 心理 学 、 管 理学 、 社 会 学 等 学 科 中 重要 
的 统计 工具 之 一 , 可 以 用 于 验证 性 因子 分 析 、 共 同方 法 偏差 检验 、 中 介 / 调 节 效 应 分 析 、 交 叉 
滞后 分 析 等 多 种 场景 。 有 研究 者 指出 ， 相 较 于 一 般 的 回归 分 析 ,， 结构 方程 模型 能 够 更 好 地 控 
制 测量 误差 , 也 支持 构建 复杂 的 多 变量 模型 ( 王 阳 等 ,2020)。 然 而 , 在 使 用 结构 方程 模型 进行 
分 析 时 ,， 却 较 少 有 研究 者 报告 统计 检验 力 的 相关 内 容 。 事 实 上 ,统计 检验 力 的 大 小 在 零 假 设 
显著 性 检验 (Null Hypothesis Significance Testing, NHST) 中 有 着 重要 的 参考 价值 。 

统计 检验 力 指 的 是 某 假设 检验 能 够 正确 地 拒绝 错误 零 假设 的 概率 。 对 统计 检验 力 进 行 分 
Br, 可 以 帮助 研究 者 确定 合理 样本 量 , 避免 造成 大 量 的 人 力 和 物力 的 浪费 ( 吴 艳 , TB, 
2011; 赵 礼 , EHF, 2019)。 此 外 ,对 评估 结构 方程 模型 拟 合 的 x? 检验 而 言 ， 零 假设 的 形式 通 
常 为 模型 拟 合 良好 或 模型 间 无 差异 , 而 备 择 假设 则 通常 为 模型 拟 合 不 良 或 模型 间 存在 差异 。 
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针对 此 类 假设 , I 类 错误 意味 着 研究 者 可 能 将 拟 合 不 展 的 模型 误 认 为 拟 合 良 好 的 模型 进行 结 
果 报 告 。 王 阳 等 人 (2020) 则 指出 , ARR BE AR, 建立 在 此 模型 基础 上 的 参数 估计 结果 都 
是 不 可 靠 的 。 在 这 种 情况 下 , I 类 错误 可 能 带 来 比 I 类 错误 更 大 的 危 

近年 来 , 有 一 些 研究 者 开始 逐渐 意识 到 使 用 结构 方程 模型 作为 统计 工具 时 统计 检验 力 
分 析 的 重要 性 , 并 开始 在 研究 中 报告 其 先 验 或 后 验 检验 力 , 或 者 至 少 在 讨论 部 分 指出 该 研究 


可 能 在 统计 检验 力 上 存 有 不 足 (Mullen & Crowe, 2017; Hollerbach et al., 2018; Kornadt et al., 
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2019; Zhang & Zaman, 2020; Zhai et al., 2021)。 但 值得 注意 的 是 ， 目 前 仍 有 许多 研究 者 在 使 用 
结构 方程 模型 时 并 未 提 及 统计 检验 力 的 有 关内 容 。 笔 者 曾 以 全 文 包含 “结构 方程 模型 ”或 
“SEM 模型 ”为 条 件 , 搜索 了 国内 心理 学 顶级 刊物 《心理 学 报 》 中 近 五 年 (2017~2021) 内 发 
表 的 论文 , 共 检 索 出 51 篇 论文 , 排除 其 中 6 篇 方法 类 研究 和 2 篇 元 分 析 研 究 , 余下 43 篇 中 
提 及 “统计 检验 力 ” 或 “统计 效力 ”的 仅 有 2 篇 研究 。 本 文 介绍 结构 方程 模型 统计 检验 力 分 
析 的 基本 原理 和 具体 分 析 流 程 ,并 进行 了 实例 演示 , 最 后 对 相关 问题 进行 了 讨论 。 
2. 结构 方程 模型 中 的 假设 检验 

假设 检验 的 技术 路 线 直 接 决 定 了 检验 力 分 析 所 能 选取 的 技术 路 线 , 结构 方程 模型 中 假 
设 检验 的 技术 路 线 大 致 可 分 为 两 种 : 基于 x 分 布 的 方法 以 及 基于 模拟 (simulation) 或 重 抽样 
(resampling) 技 术 的 方法 。 
2.1 基于 入 分 布 的 结构 方程 模型 检验 
基于 92 分 布 的 结构 方程 模型 检验 主要 包括 用 于 评价 模型 拟 合 (或 悉 套 模型 比较 ) 的 92 t 
验 、A 态 检验、 等 效 性 检验 、 基 于 模型 拟 合 指数 的 检验 以 及 用 于 评价 模型 参数 (或 参数 的 函数 ， 
如 中 介 效 应 ) 的 Wald 检验 、 拉 格 朗 日 乘 子 检验 (Lagrange Multiplier test, LM 检验 ) 和 似 然 比 检 


验 (Satorra, 1989; Chou & Bentler, 1990; Gonzalez & Griffin, 2001; Yuan & Chan, 2016; Yuan et 
al., 2016; 王 阳 等 ,2020)。 这 些 检 验 涵盖 了 结构 方程 模型 相关 假设 检验 的 方方面面 甚至 可 以 
说 , 目前 与 结构 方程 模型 有 关 的 参数 检验 方法 要 么 本 身 就 是 x 检验 , 要 么 可 以 通过 适当 变形 
转化 为 凶 检 验 。 根 据 构建 统计 量 时 的 细微 差别 ， 上述 检 验 又 可 被 划分 为 三 个 子 类 。 

首先 是 闻 检 验 、AX? 检验、 等 效 性 检验 和 似 然 比 检验 , 这 四 种 检验 都 以 结构 方程 模型 的 
拟 合 函数 为 基础 。 结 构 方程 模型 的 拟 合 函数 衡量 了 样本 协 方差 矩阵 S 与 假设 模型 对 应 的 协 
方差 矩阵 Zo 之 间 的 差距 , 亦 即 变 量 间 关 系 在 假设 模型 与 实际 数据 中 的 差距 ( 吴 明 隆 ,2010; E 
BHSE, 2020)。Muthén(2004) 指 出 ， 当 样本 来 自 假设 模型 所 对 应 的 正 态 总 体 时 , 拟 合 函 数值 的 
N-1 倍 (N 为 样本 量 ) 随 着 N 的 增 大 渐 近 服从 中 心 化 闻 分 布 。 另 有 学 者 则 指出 ， 当 假定 样本 来 
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自 与 假设 模型 差距 小 于 给 定 值 ( 记 作 Fo) 的 正 态 总 体 时 ， 此 时 拟 合 函 数值 的 N-7 倍 渐进 服从 


非 中心 参 数 为 (V-7)Fo 的 非 


Chan, 2016)。 据 此 ， 可 以 构建 出 零 假 设 成 立时 ?检验 和 等 


Ay? 检验 和 似 然 比 检验 则 是 在 前 者 的 大 
约束 条 件 前 后 实际 数据 与 假设 模型 的 差距 是 否 显著 增 大 来 判定 添加 的 约束 条 


件 是 否 “ 合 理 ”。 


这 类 检验 与 第 一 类 检验 差别 不 大 ,其 
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二 类 是 基于 模型 拟 合 指 数 的 检验 。 
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E, 增加 了 对 假设 模型 的 约束 条 伯 


效 性 检验 所 对 应 的 统计 量 


从 原理 上 讲 ， 


Hut» 和 2 分布 (Satorra & Saris, 1985; MacCallum et al., 1996; Yuan & 


分 布 。 而 


F, 通过 比较 添加 


本 质 是 通过 拟 合 指数 的 定义 式 将 拟 合 指数 换算 为 非 中 心 参数 , 因此 该 方法 仅 能 用 于 根据 非 


中 心 参数 定义 的 拟 合 指数 (如 RMSEA、CFI、MFI、GFL 等 )。 此 类 检验 有 两 种 相似 但 略 有 不 


同 的 技术 路 线 ， 一 种 将 基于 模型 拟 合 指数 的 零 假 设 换算 为 基于 非 中 心 参数 的 零 假 设 , 选择 非 
中 心 参数 的 临界 值 作为 构建 ?统计 量 分 布 的 依据 , MacCallum 等 人 (1996) 和 MacCallum 等 人 


(2006) 就 采用 了 这 种 做 法 。 另 一 种 则 根据 Venables(1975) 的 研究 ， 以 样 
估计 值 为 依据 构建 总 体 非 


而 得 到 拟 合 指数 的 1-a 置信 区 间 ， à 


的 等 效 性 检验 采用 了 此 类 做 法 (Yuan et al., 2016; Marcoulides & Yuan, 2017; Ef 


第 三 类 是 Wald 检验 和 LM 检验 。 


计 , 构建 z 统计 量 进行 检验 。 
数 的 估计 和 检验 。 当 处 理 参数 函数 形式 的 假设 时 , Wald 检验 对 参数 函数 在 其 真 值 邻 域内 进行 


一 阶 Taylor 展开 求 得 其 线性 近似 值 ， 进 而 完成 对 待 检验 指标 的 标准 误 估计 (Chou & Bentler, 
1990; Gonzalez & Griffin, 2001)。 此 时 Wald 检验 亦 称 为 Delta Yk, 叶 宝 娟 和 温 


过 如 何 利 用 该 方法 对 多 维 
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以 零 假 设 成 立时 的 模型 为 


心 参数 的 1-a 置信 区 间 ， 


测验 合 


基础 


， 将 假设 条 件 重 写 为 约束 函 


原 假设 检验 问题 变 为 检验 也 


乘 数 为 0， 则 ; 


再 根据 定义 式 将 其 换算 
而 完成 统计 推断 。Yuan 等 人 提出 的 基于 CFI 与 RMSEA 


成 信 度 的 置信 区 间 进 


2.2 基于 模拟 或 重 抽样 技术 的 结构 方程 模型 检验 


近年 来 ， 随 着 计算 机 技术 的 发 展 ,模拟 技术 和 重 抽样 
这 些 技术 被 广泛 地 运用 在 包括 参数 佑 计 与 假设 检验 在 内 
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意味 着 新 假设 条 件 的 加 入 没有 改变 原 假设 模型 。 


忠 麟 (2012) 曾 介 


本 计算 出 的 非 中 心 参数 
回 拟 合 指数 ， 从 


等 , 2020)。 
与 前 两 类 方法 不 同 的 是 ， 此 类 方法 不 再 关注 样 
合 函 数值 本 身 ( 非 中 心 参 数 也 是 关于 拟 合 函 


本 的 拟 


数 的 函数 ), 而 是 通过 对 关键 参数 及 其 标准 误 的 估 


其 中 Wald 检验 以 备 择 假设 成 立时 的 模型 为 基础 完成 对 目标 参 


行 估计 。 与 Wald 检验 相对 ,LM 检验 
数 的 形式 ， 并 构建 拉 格 朗 日 


函数 将 


立 格 朗 日 乘 数 是 否 为 0 的 问题 (Chou & Bentler, 1990)。 若 拉 格 朗 日 


技术 越 来 越 多 地 出 现在 统计 学 领域 ， 
场合 .不 同 于 基于 理论 分 布 的 


2 本 文中 出 现 的 “GFI” 和 “AGFI” 指 的 分 别 是 “Unbiased goodness of fit index” 和 “Unbiased adjusted 


goodness of fit index” 


, IE “Goodness of fit index” FI “Adjusted goodness of fit index” 


， 详 见 补充 材料 。 


chinaXiv:202110.00074v2 


经 典 方法 ， 此 类 方法 通常 以 模拟 或 重 抽样 的 方式 得 到 目标 统计 量 的 经 验 分 布 , 并 以 经 验 分 布 


代 蔡 理论 分 布 进行 统计 推断 ， 


因此 往往 用 于 目标 统计 量 的 理论 分 布 不 明确 或 难于 计算 的 情 


形 。 此 类 方法 中 最 为 重要 的 两 种 分 别 是 基于 模拟 技术 的 Monte Carlo (MC) 方 法 3 与 基于 重 抽 


样 技术 的 Bootstrap 方法 。 


MC 方法 最 早 由 Stanistaw Ulam 等 人 于 20 世纪 40 年 代 中 期 提出 ， 其 核心 ， 
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时 想 为 采用 模 


EL” ” 某 一 概率 模型 所 描绘 的 过 程 ， 以 模拟 得 到 的 结果 作为 目标 问题 的 近似 解 。 


MacKinnon 等 人 (2004) 在 评估 中 介 效 应 检验 的 M 法 和 经 验 M 法 时 同时 引入 了 MC 法 。 


Preacher #0 S 


elig(2012) 进 一 步 指 出 MC 法 不 仅 可 以 用 于 线性 回归 和 路 径 分 析 模 


型 ， 也 可 直接 


扩展 至 潜 变量 模型 。 该 方法 假定 模型 参数 服从 联合 正 态 分 布 ， 其 均值 向 量 与 协 方差 矩阵 通常 


过 程 ， 以 此 为 基础 


夫 链 蒙特 卡 洛 (Markov chain Monte Carlo, MCMC) 方 法 ， 此 类 方法 主要 月 


由 对 原始 模型 的 极 大 似 然 估计 确定 。 MC 法 将 模型 参数 视 为 从 确 
进行 随机 抽样 ,最 终 基 于 抽样 结果 构建 中 介 效 应 的 经 验 分 布 并 进行 统计 推 
Wt(Preacher & Selig, 2012)。 此 外 ,在 结构 方程 建 模 中 还 有 一 类 Monte Carlo 方法 一 一 马尔 可 
于 贝 叶 斯 结构 方程 


定 的 联合 正 态 分 布 中 抽样 的 


建 模 中 。 由 于 对 “参数 ”这 一 重要 概念 的 基础 假设 不 同 ， 贝 叶 斯 统计 中 假设 检验 和 统计 检验 
力 分 析 的 过 程 均 与 基于 频率 学 派 的 传统 统计 学 有 着 较 大 的 区 别 , 本文 的 内 容 仍然 在 频率 学 


派 的 框架 内 ， 


故 对 该 方法 不 做 过 多 介绍 。 


注意 到 ， 基 于 模拟 技术 的 结构 方程 模型 检验 仍然 对 分 布 有 着 较 强 的 假设 , MC 法 假定 路 
径 系数 服从 联合 正 态 分 布 ,基于 贝 叶 斯 统计 的 MCMC 法 需要 给 出 模型 参数 的 先 验 分 布 与 给 
定 参 数 时 样本 的 条 件 分 布 。 相 比 之 下 ,基于 重 抽样 技术 的 检验 对 样本 分 布 的 假设 更 弱 ， 
Efron(1979) 提 出 的 Bootstrap 法 仅 要 求 样本 是 从 某 一 确定 总 体 中 随机 抽样 得 到 
形态 并 无 更 多 要 求 。 不 仅 如 此 , Bootstrap 法 同样 基于 “利用 经 验 分 布 代 蔡 理 论 分 布 进行 统计 


推断 ”的 思路 , 在 面 对 一 些 理论 分 布 不 明 或 难以 计算 的 场合 表现 良好 。 由 于 Bootstrap 法 具 


有 前 述 技术 优势 ,近年 来 ,该 方法 已 逐渐 成 为 结构 方程 模型 假设 检验 问题 中 的 “ 显 学 ”, 被 


的 ， 对 其 分 布 


广泛 应 用 于 多 种 不 同 场 合 。Zhang 和 Savalei(2016) 指 出 可 以 用 Bootstrap 法 估计 结构 方程 模 
型 拟 合 指数 的 置信 区 间 。Rosseel 在 其 编写 的 及 程序 包 lavaan 中 给 出 了 求 结构 方程 模型 路 径 
系数 的 Bootstrap 标准 误 的 程序 (Rosseel, 2021)。 更 有 研究 者 利用 Bootstrap 法 估计 一 些 根据 


3 Monte Carlo 方法 实际 上 代表 的 是 以 计算 机 模拟 为 核心 的 一 组 方法 , 为 了 和 后 文 统 计 检 验 力 分 析 的 Monte 


Carlo 方法 加 以 


区 分 , 假设 检验 的 Monte Carlo 方法 本 文 统一 简称 为 MC 方法 /MC 法 , 特此 说 明 。 


4 事实 上 , Bootstrap 法 还 分 为 参数 Bootstrap 法 和 非 参数 Bootstrap 法 ,参数 Bootstrap 法 由 于 需要 了 解 随机 
变量 的 分 布 型 态 , 在 实际 中 使 用 不 多 。 本 文中 若 不 加 说 明 , Bootstrap 法 特 指 非 参数 Bootstrap 法 。 


结构 方程 模型 路 径 系数 计 外 


得 到 的 更 复杂 的 指标 。 屠 金 路 等 人 (2005) 利 用 Bootstrap 法 对 问 


卷 合 成 信 度 的 置信 区 间 进 行 了 估计 ; 张 涵 和 康 飞 (2016) 介 绍 了 基于 Bootstrap 的 多 重 中 介 效 


应 分 析 。 方 杰 等 人 (2011) 的 研究 及 温 忠 鹿 和 叶 宝 娟 (2014a, 2014b) 的 研究 则 分 别 介绍 了 


Bootstrap 法 在 中 介 效 应 与 有 调节 的 中 介 效 应 检验 中 的 应 用 。 


表 1 中 按照 不 同类 别 总 结 了 结构 方程 模型 涉及 到 的 常见 检验 的 应 用 场景 、 零 假设 、 统 计 


量 构建 方式 、 零 假设 成 立时 统计 量 对 应 的 分 布 以 及 使 用 中 的 优 缺 点 。 


注 : 由 于 加 入 内 容 后 表格 过 大 ,难以 在 本 页 显示 ,我 们 在 修改 时 将 更 新 后 的 表 1 置 于 文 末 


3， 结 构 方程 模型 假设 检验 统计 检验 力 的 分 析 方 法 

统计 检验 力 的 分 析 方法 与 分 析 对 象 所 采用 的 技术 路 径 息息相关 ,结构 方程 模型 假设 检 
验 统计 检验 力 的 分 析 方法 也 可 大 致 分 为 两 条 不 同 的 技术 路 线 : 基于 9? 分 布 的 方法 和 基于 模 
拟 的 方法 。 其 中 , 基于 允 分 布 的 方法 可 以 用 于 解决 使 用 经 典 方法 进行 结构 方程 模型 假设 检验 


时 的 统计 检验 力 分 析 问 题 ， 
的 统计 检验 力 分 析 问 题 。 


而 基于 模拟 的 方法 可 以 用 于 解决 一 般 的 结构 方程 模型 假设 检验 


3.1 基于 多 分 布 的 统计 检验 力 分 析 
顾名思义 ,基于 入 分 布 的 方法 以 巡 分 布 的 相关 理论 为 基础 ,此 类 方法 可 以 追溯 到 Satorra 


和 Saris 早年 的 一 系列 研究 (Satorra & Saris, 1985; Satorra, 1989; Satorra et al., 1991). 


Satorra(1989) 指 出 ， 当 前 述 假设 检验 的 加 不 成 立时 ， 其 对 应 统计 量 (对 Wald 检验 和 LM 检验 


而 言 是 统计 量 的 平方 ) 服 从 非 中 心 e 分布 且 其 相应 的 非 中 心 参数 可 以 通过 总 体 协 方差 矩阵 
真 值 拟 合 零 假设 模型 计算 出 的 统计 量 得 到 。 进 而 可 以 通过 计算 备 择 假设 成 立 条 件 下 目标 统计 
落 入 拒绝 域 的 概率 求 得 给 定 检验 的 统计 检验 力 。 设 零 假 设 和 备 择 假设 成 立时 对 应 的 非 中 心 


te 


参数 分 别 为 和 51， 则 统计 检验 力 1-B 由 下 式 给 出 : 


1 一 B = P{X3.s, < Twzc]( 拒 绝 域 在 左 侧 , 对 应 等 效 性 检验 ) (1a) 


或 1 一 B = P{x376, > Twic}( 拒 绝 域 在 右 侧 ,对 


El 


除 等 效 性 检验 外 的 其 他 检验 ) (1b) 


其 中 x&ps 表 示 自 由 度 为 df 非 中 心 参数 为 6 的 非 中 心 卡 方 分 布 , Tuc 由 下 式 确定 : 


p{X34s, > Tuc} = a( 拒 绝 域 在 左 侧 ,对 应 等 效 性 检验 ) (2a) 


或 Pfx346, < Ture} = a( 拒 绝 域 在 右 侧 ,对 应 除 等 效 性 检验 外 的 其 他 检验 ) (28) 


由 于 3 在 Ho HIERN CARA ASE, PALE, 基于 入 分 布 的 统计 检验 力 分 析 的 关键 在 于 估计 61。 
根据 估计 非 中 心 参数 的 方法 不 同 , 此 类 方法 又 可 分 为 利用 总 体 协 方差 矩阵 三 进行 估计 的 


Satorra-Saris 法 和 利用 拟 合 指 数 ( 如 RMSEA、CFI、MFI 等 ) 进 行 估计 的 MacCallum 法 。 


3.1.1 Satorra—Saris 法 


N 
if 
I 


Satorra—Saris 法 主要 基于 Satorra 和 Saris 先前 的 研究 ， 


其 具体 步骤 如 下 (Satorra & Saris, 


1985; Satorra et al., 1991): 
第 一 步 : 定义 备 择 假设 成 立 情 况 下 的 模型 ， 指 定 备 择 假设 成 立 情况 下 的 理论 总 体 协 方差 
FEM X. 
第 二 步 : 设 定 样本 容量 为 真实 样本 容量 N, 利用 三 拟 合 零 假设 模型 。 
第 三 步 : 将 第 二 步 中 计算 得 到 的 统计 量 (针对 x 分 布 的 统计 量 为 统计 量 本 身 , 针对 标准 
正 态 分 布 的 统计 量 为 统计 量 的 平方 ) 作 为 61 的 估计 , 代入 (1a) 或 (lb) 完 成 统计 检验 力 计算 。 
使 用 Satorra-Saris 法 估计 达到 某 统 计 检 验 力 所 需 样 本 量 的 过 程 与 估计 统计 检验 力 类 似 ， 
只 是 第 二 步 中 需要 反复 设 定 样本 容量 N 进行 多 次 估计 ,并 计算 不 同样 本 容量 对 应 的 统计 检 
验 力 , 最 终 从 中 选取 满足 目标 统计 检验 力 所 需 的 最 小 样本 容量 

Satorra-Saris 法 通过 定义 不 同 于 零 假 设 模 型 的 备 择 模型 重 构 总 体 协 方差 矩阵 ， 并 据 此 计 
算 统计 检验 力 。 其 中 , 对 备 择 模型 进行 定义 的 过 程 不 仅 需 要 研究 者 给 出 明确 的 模型 结构 ， 还 
需要 研究 者 指定 大 量 的 参数 (如 路 径 系 数 、 测 量 残 差 等 )。 因 此 ， 从 应 用 的 角度 考虑 , Satorra— 
Saris 法 更 适合 在 研究 者 对 备 择 模 型 的 了 解 程度 较 深 的 情况 下 使 用 (如 针对 模型 参数 的 检验 )。 


3.1.2 MacCallum 法 


Satorra-Saris 法 对 备 择 模型 过 强 的 假设 从 客观 上 提高 了 该 方法 的 使 用 难度 ,也 限制 了 其 


适用 范围 。MacCallum、Browne 和 Sugawara(1996) 给 出 了 另外 一 种 估计 非 中 心 参数 的 方法 ， 


该 方法 同样 基于 Satorra 和 Saris(1985) 给 出 的 公式 , 但 在 估计 Furo BT, MacCallum 等 人 给 出 了 
一 个 基于 RMSEA 的 估计 方法 。 他 们 指出 ,使 用 样本 拟 合 函数 估计 值 Fu 作为 对 总 体 参数 
Furo 的 估计 是 有 偏 的 ， 且 Fwo 与 Fu 之 间 有 如 下 关系 (MacCallum et al., 1996): 


Furo = Fi a g - RMSEA? (3) 
MLO ~ "ML yy f 
据 此 , MacCallum 等 人 提出 基于 RMSEA 指标 的 统计 检验 力 分 析 步 又 (MacCallum et al., 1996): 
第 一 步 : 根据 零 假 设 和 备 择 假设 成 立时 对 应 的 RMSEA 值 、 模型 自由 度 以 及 样本 量 分 别 
计算 出 零 假 设 和 备 择 假设 成 立时 的 非 中 心 参数 50 和 61。 


第 二 步 : 代入 (1a) 或 (lb) 完 成 统计 检验 力 计 算 。 


D 


与 Satorra-Saris 法 类 似 ， 当 使 用 MacCallum 法 估计 达到 某 统计 检验 力 所 需 样本 量 时 ,也 
只 需要 在 第 一 步 中 反复 设 定 不 同 的 N 对 统计 检验 力 进行 多 次 估计 , 最 终 从 中 选取 满足 目标 
统计 检验 力 所 需 的 最 小 样本 容量 。 

Kim(2005) 对 MacCallum 法 进行 了 进一步 简化 和 扩展 , 给 出 了 一 个 快速 求 取 给 定 统 计 检 
验 力 情况 下 所 需 样本 量 的 公式 


6 
N 二 一 
RMSEA = Gf RMSEA? + 


其 中 df KIW A hE, RMSEA 为 备 择 假设 成 立时 对 应 的 RMSEA 值 , 8 为 显著 性 水 平 为 a、 
统计 检验 力 为 1-B 情况 下 对 应 的 非 中 心 参数 。 非 中 心 参数 可 以 通过 查 表 或 者 使 用 统计 软件 
计算 得 到 , Kim(2005) 在 他 们 的 文章 最 后 同时 给 出 了 Ho: Furo = 0 时 利用 计算 机 快速 求 取 5 的 
SPSS 与 SAS 脚本 。 此 外 , Kim(2005) 在 他 的 文章 中 还 提出 了 基于 CFI、MFI 与 GFI 的 计算 公 
式 。 

相 比 于 Satorra-Saris 法 , MacCallum 法 采用 拟 合 指数 对 模型 的 拟 合 情 况 进行 了 某 种 意义 
上 的 “打包 ”, 即 利用 拟 合 指数 对 假设 模型 与 实际 数据 的 差距 进行 了 概括 , 使 研究 者 免 于 设 
定 复杂 的 备 择 模型 ,在 分 析 模 型 拟 合 检验 的 统计 检验 力 时 , MacCallum 法 仅 需 要 模型 自由 度 、 
零 假设 和 备 择 假设 成 立时 对 应 的 拟 合 指数 值 [ 若 对 关于 CFI 的 检验 进行 检验 力 分 析 时 还 需要 
提供 基础 模型 (base model) 的 拟 合 函 数值 mas 与 自由 度 dfs, 关于 GFI 的 检验 还 需要 提供 协 
方差 矩阵 的 秩 刀 ,大 大 简化 了 统计 检验 力 的 分 析 难 度 。 但 由 于 模型 参数 与 拟 合 指数 之 间 存 在 
杂 的 非 线性 关系 , MacCallum 法 无 法 分 析 模 型 参数 检验 的 统计 检验 力 。 
3. 2 基于 模拟 的 统计 检验 力 分 析 

采用 模拟 方法 对 结构 方程 模型 假设 检验 进行 统计 检验 力 分 析 早 年 可 追溯 到 Muthén 与 
Muthén(2002) 的 研究 。 采用 模拟 方法 进行 统计 检验 力 分 析 的 思路 与 采用 模拟 方法 进行 结构 方 
程 模型 假设 检验 的 思路 相似 , 即 利用 模拟 的 方式 确定 P( 拒 绝 Ho19 = 0 ) 的 经 验 值 。 具体 而 言 ， 
此 类 方法 通常 采用 模拟 的 方式 从 参数 为 9" 的 特定 分 布 中 生成 多 组 随机 样本 , 而 后 针对 每 组 
随机 样本 执行 同样 的 假设 检验 程序 , 最 后 统计 全 部 已 进行 的 假设 检验 中 拒绝 零 假设 的 频率 ， 
以 此 作为 对 该 检验 统计 检验 力 的 分 析 结 果 。Muthén 与 Muthén(2002) 的 研究 给 出 了 利用 Monte 
Carlo 模拟 计算 基于 经 典 方法 的 结构 方程 模型 假设 检验 的 统计 检验 力 的 分 析 过 程 与 Mplus 程 
序 。Zhang(2014) 的 研究 指出 Monte Carlo 模拟 法 同样 可 以 估计 基于 Bootstrap 方法 的 中 介 效 
应 检验 的 统计 检验 力 , 并 给 出 了 一 个 基于 R 语言 bmem 程序 包 的 示例 程序 。 Schoemann 等 人 
(2017) 的 研究 在 实现 对 中 介 模 型 的 统计 检验 力 分 析 时 使 用 了 Monte Carlo 模拟 法 估计 基于 
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MC 法 的 中 介 效 应 检验 的 统计 检验 力 。 总 的 来 说 , 使 用 Monte Carlo 模拟 法 进行 统计 检验 力 
分 析 的 步骤 大 体 如 下 (Muthkn & Muthén, 2002; Zhang, 2014; Schoemann et al., 2017): 

第 一 步 : 根据 备 择 假 设 成 立时 的 情况 设 定 用 于 生成 模拟 数据 的 均值 向 量 (模型 不 考虑 均 
值 结构 时 为 零 向 量 ) 与 协 方差 矩阵 。 此 步骤 中 既 可 以 如 Schoeman 等 人 (2017) 直 接 指 定 协 方 
差 矩阵 ， 也 可 如 Zhang(2014) 通 过 定义 模型 来 间接 指定 协 方差 矩阵 。 此 外 , 如 果 研 究 者 希望 模 
拟 非 正 态 分 布 的 数据 ,通常 还 需要 给 出 目标 分 布 的 峰 度 和 偏 度 。 


第 二 步 : 确定 需要 进行 Monte Carlo 模拟 的 重复 轮 数 R 与 样本 容量 W。Zhang(2014) 建 议 


R 通常 不 应 少 于 1000。 

第 三 步 : 在 每 一 轮 Monte Carlo 模拟 中 , 重复 执行 下 述 步骤 : (1) 根 据 给 定 的 分 布 参数 与 
样本 量 生成 服从 特定 分 布 的 随机 数 作为 模拟 样本 ; (2) 对 该 模拟 样本 进行 结构 方程 建 模 并 进 
行 检 验 ， 有 具体 的 检验 方法 可 以 是 表 1 中 的 任意 一 种 。 

第 四 步 : 统计 出 每 一 轮 统计 检验 的 结果 , 设 在 R 轮 检验 中 有 /7 轮 结果 拒绝 了 原 假 设 ， 则 


最 终 的 统计 检验 力 为 


此 外 , 为 了 说 明 Monte Carlo 模拟 结果 的 准确 性 , Muthén 和 Muthén(2002) 建 议 研究 者 在 
报告 结果 的 同时 报告 结构 方程 模型 中 各 参数 的 理论 真 值 、 参 数 估 计 结 果 的 平均 值 和 标准 差 、 
参数 估计 标准 误 的 平均 值 以 及 参数 95% 置 信 区 间 对 理论 真 值 的 覆盖 率 。Muthén 和 
Muthén(2002) 进 一 步 提出 衡量 Monte Carlo 模拟 结果 是 否 可 靠 的 标准 : (1) 模 型 中 任何 参数 及 
其 标准 误 真 值 与 估计 值 之 间 的 相对 误差 不 超过 10% (参数 估计 结果 的 标准 差 视 为 标准 误 的 真 
值 ); (2) 模 型 中 关键 参数 及 其 标准 误 真 值 与 估计 值 之 间 的 相对 误差 不 超过 5%; (3) 参 数 的 95% 
© 置信 区 间 对 理论 真 值 的 覆盖 率 应 当 在 0.91 至 0.98 之 间 。 

当 使 用 Monte Carlo 法 估计 达到 某 统计 检验 力 所 需 样本 量 时 ， 也 只 需要 在 第 一 步 中 反复 

设 定 不 同 的 W 对 统计 检验 力 进行 多 次 估计 ,而 后 通过 线性 插值 法 对 所 需 样 本 量 进行 估计 。 与 

基于 % 分 布 的 统计 检验 力 分 析 方 法 相 比 , 结构 方程 模型 统计 检验 力 的 Monte Carlo 分 析 由 于 

不 需要 计算 目标 统计 量 的 理论 分 布 , 因此 可 以 广泛 适用 于 对 各 种 检验 进行 统计 检验 力 分 析 

的 情形 , 在 给 出 峰 度 和 偏 度 时 ， 该 方法 还 可 以 在 一 定 程度 上 分 析 总 体 分 布 为 非 正 态 分 布 的 情 

形 。 但 结构 方程 模型 统计 检验 力 分 析 的 Monte Carlo 法 需要 研究 者 至 少 给 出 用 以 生成 模拟 数 

据 的 协 方差 矩阵 ， 这 意味 着 研究 者 需要 对 备 择 假 设 成 立 下 数据 的 分 布 情况 有 足够 的 认识 。 另 

外 ,该 方法 本 身 对 研究 者 的 编程 水 平 要 求 较 高 ,， 且 计算 量 较 大 , 计算 起 来 较为 耗 时 , 在 一 定 
程度 上 影响 了 该 方法 的 使 用 和 推广 。 


结构 方程 模型 统计 检验 力 分 析 的 程序 资源 

为 了 方便 研究 者 对 统计 检验 力 进行 分 析 , 多 名 研究 者 提供 了 在 线 或 离线 的 相关 程序 。 基 
于 和 2 分 布 的 统计 检验 力 分 析 方 法 中 , Preacher 和 Coffman(2006) 提 供 了 一 个 基于 MacCallum 
去 的 统计 检验 力 在 线 计 算 程 序 , 可 以 计算 给 定 样本 量 情况 下 ,单个 模型 拟 合 或 庶 套 模型 拟 合 
差异 的 统计 效力 或 达到 某 统计 效力 所 需 的 样本 量 .Moshagen 提供 了 一 个 RR 程序 包 semPower, 
该 程序 包 基 于 Moshagen 和 Erdfelder(2016) 的 研究 ,可 以 计算 基于 Furos RMSEA, McDonald 
拟 合 指数 、GFI 与 AGFI 的 统计 检验 力 ， 当 指定 零 假设 和 备 择 假设 成 立 情况 下 对 应 的 协 方差 
矩阵 时 ， 该 程序 包 也 可 使 用 Satorra-Saris 法 分 析 统 计 检 验 力 。Jak 等 人 (2021) 给 出 了 一 个 基于 
Shiny 开发 的 在 线程 序 power4SEM( 网 址 为 https://sjak.shinyapps.io/power4SEM/), 该 程序 同时 
LİF Satorra-Saris 法 和 MacCallum 法 , Jak 等 人 在 他 们 的 论文 中 详细 介绍 了 该 程序 的 使 用 方 
式 。 
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基于 模拟 的 统计 检验 力 分 析 方法 中 , Muthén 和 Muthén(2002) 提 供 了 使 用 Monte Carlo 法 
分 析 结 构 方程 模型 统计 检验 力 的 Mplus 程序 框架 。Zhang 和 Wang(2013) 提 供 了 可 以 用 来 分 
析 中 介 效 应 统计 检验 力 的 及 语言 程序 包 bmem。 该 程序 包 使 用 Monte Carlo 法 进行 统计 检验 
力 分 析 ， 可 以 对 中 介 效 应 的 Sobel 检验 与 Bootstrap 检验 的 统计 检验 力 进行 分 析 。 特 别 的 是 ， 
该 程序 包 还 提供 给 研究 者 自主 设 定 分 布 的 峰 度 与 偏 度 的 功能 以 模拟 总 体 为 非 正 态 分 布 时 的 
情形 。Schoemann 等 人 (2017) 给 出 了 一 个 使 用 Monte Carlo 法 分 析 中 介 效 应 统计 检验 力 的 在 
线程 序 (网 址 为 https:/schoemanna.shinyapps.io/mc power med/)， 该 程序 基于 Shiny FR, 可 
用 来 处 理 简单 中 介 模 型 、 两 个 中 介 变 量 的 平行 中 介 模 型 和 两 个 中 介 变 量 的 链 式 中 介 模 型 三 种 
情况 下 的 中 介 效 应 统计 检验 力 分 析 。Wang 和 Rhemtulla(2021) 则 给 出 了 一 个 分 析 结 构 方 程 模 
型 参数 检验 统计 检验 力 的 在 线程 序 (网 址 为 : https://yilinandrewang.shinyapps.io/pwrSEM/)。 

表 2 中 总 结 了 结构 方程 模型 统计 检验 力 分 析 方 法 的 适用 场合 、 需 要 的 计算 量 以 及 相应 
的 程序 资源 。 


2 结构 方程 模型 统计 检验 力 分 析 方 法 
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4， 应 用 实例 

下 面 用 一 个 例子 演示 在 实际 研究 中 如 何 灵 活 使 用 各 种 方法 对 结构 方程 模型 的 统计 检验 
力 进 行 分 析 。 注 意 到 ,统计 检验 力 的 分 析 包 括 先 验 分 析 和 后 验 分 析 ,， 二 者 对 应 研究 过 程 中 的 
不 同时 期 ， 且 不 能 相互 替代 ( 吴 艳 ,， 温 忠 麟 ,2011; 赵 礼 ， 王 晖 ,2019)。 但 先 验 分 析 与 后 验 分 析 
在 分 析 步 骤 上 差异 不 大 , 且 亦 有 研究 者 指出 , 在 NHST 中 , 统计 检验 力 是 值 的 函数 , 在 已 
知 p 值 的 情况 下 统计 检验 力 没有 给 出 额外 信息 ,因此 建议 无 需 进行 后 验 分 析 (Hoenig & 
Heisey, 2001)。 考 虑 到 篇 幅 所 限 ， 本 文 仅 展 示 先 验 分 析 的 过 程 ， 后 验 分 析 的 部 分 参见 补充 材 
料 。 演 示 的 样 例 来 自 Zhai 等 人 (2021) 的 研究 , 该 研究 调查 了 新 冠 疫 情 期 间 个 体 的 心理 健康 与 
创伤 后 成 长 的 基本 情况 ， 并 讨论 了 在 控制 疫情 相关 压力 和 年 龄 的 情况 下 , 情绪 创造 力 对 个 体 
心理 健康 的 保护 作用 以 及 对 个 体 创伤 后 成 长 的 促进 作用 和 可 能 的 作用 机 制 . 具 体 的 模型 设置 
参见 图 1, 统计 软件 采用 R 3.6.3 的 lavaan 程序 包 和 semPower 程序 包 , 本 部 分 全 部 R 代码 参 
见 补充 材料 部 分 。 


压力 后 成 长 
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图 1 样 例 的 结构 方程 模型 示意 图 (虚线 表示 该 变量 为 控制 变量 ) 

从 检验 目的 的 角度 出 发 , 使 用 结构 方程 模型 进行 验证 性 因子 分 析 的 主要 目的 是 通过 对 
模型 拟 合 的 评价 考察 样 本 数据 是 否 符 合 量 表 所 的 理论 建构 ， 即 考察 量 表 题 项 与 测量 因子 间 
的 关系 是 否 同 研究 者 的 假设 所 一 致 ; 而 使 用 结构 方程 模型 进行 中 介 效 应 分 析 , 特别 是 对 显 变 
量 链 式 中 介 模 型 进行 中 介 效 应 分 析 时 ,研究 者 的 目的 则 更 多 指向 对 路 径 系 数 函 数 ( 即 中 介 效 
应 ) 的 分 析 与 检验 。 此 外 , 注意 到 图 1 中 的 模型 为 饱和 模型 ， 有 研究 者 指出 ， 此 时 模型 为 完美 
拟 合 ， 可 不 关注 模型 拟 合 指数 ,， 仅 讨论 路 径 系数 即 可 ( 张 莉 等 ,2019)。 综 合 考虑 检验 目的 与 实 


际 情况 ,对 该 研究 的 统计 检验 力 进 行 分 析 时 , 需要 对 验证 性 因子 分 析 部 分 模型 拟 合 检验 的 统 
计 检 验 力 与 链 式 中 介 模 型 部 分 模型 参数 检验 的 统计 检验 力 进行 分 析 , 链 式 中 介 效 应 部 分 选 
择 “ 情 绪 创 造 力 一 领悟 社会 支持 一 情绪 调节 自我 效能 感 一 压力 后 成 长 ”路 径 的 中 介 效 应 值 为 
分 析 对 象 。 
4.1 模型 拟 合 检验 的 先 验 检验 力 分 析 
4.1.1 使 用 Satorra-Saris 法 分 析 模 型 拟 合 检验 的 先 验 检验 力 

在 进行 先 验 检 验 力 分 析 时 ,首先 要 设 定 研 究 对 象 理 论 上 的 总 体 协 方差 矩阵 (下 简称 理论 
协 方差 矩阵 )， 即 备 择 假设 成 立时 对 应 的 协 方差 矩阵 。 与 上 检验 、 方 差分 析 的 先 验 检 验 力 分 析 
类 似 , 在 结构 方程 模型 相关 检验 的 先 验 检 验 力 分 析 中 , 理论 协 方差 矩阵 也 可 以 通过 参考 前 期 
理论 、 己 发 表 的 类 似 成 果 或 研究 者 提出 的 标准 所 确定 ( 吴 艳 , 温 忠 鹿 ,2011)。 通 常情 况 下 ， 建 
议 研究 者 优先 参考 已 有 类 似 研 究 中 报告 的 结果 设 定理 论 协 方差 矩阵 ， 当 缺乏 既往 研究 或 既 
往 研究 信息 不 充足 时 , 研究 者 可 以 参考 前 期 理论 或 其 他 研究 者 提出 的 标准 所 确定 。 此 外 , Jak 
等 人 (2021) 建 议 研究 者 在 设 定理 论 协 方差 矩阵 时 ， 可 以 以 标准 化 的 形式 指定 参数 值 ， 并 且 参 
考 相关 研究 领域 中 “小 、 中 、 大 ”所 对 应 的 经 验 值 代入 待定 参数 。 同 时 ， 由 于 通常 情况 下 确 
定 协 方差 矩阵 所 需要 的 参数 数量 较为 庞大 , Jak 等 人 (2021) 还 建议 研究 者 在 设 定理 论 协 方差 
和 矩阵 时 需要 考虑 那些 在 理论 上 能 够 影响 拒绝 或 接受 Hy 的 重要 参数 , 例如 路 径 系 数 、 误 差 相 
关 、 潜 变量 相关 等 。 
通过 查询 资料 发 现 : 既往 研究 中 , 文书 峰 等 人 (2009) 的 研究 较为 全 面 地 报告 了 可 以 用 于 
生成 理论 协 方差 矩阵 的 全 部 参数 ,， 且 该 研究 样本 量 较 大 (验证 性 因子 分 析 部 分 N=818)。 因 此 ,， 
在 之 后 的 分 析 中 ， 以 该 模型 为 基准 生成 理论 协 方差 矩阵 。 本 例 中 控制 I 类 错误 水 平 a 为 0.05， 
求 得 统计 检验 力 为 0.95 时 对 应 的 样本 量 约 为 134 人 。 此 外 ,对 不 熟悉 R 的 研究 者 ， 也 可 以 使 
常用 的 结构 方程 建 模 软件 (如 LISREL、EQS、AMOS、Mplus 等 ) 结 合 G*Power 进行 分 析 ， 
具体 的 步 又 参见 补充 材料 部 分 。 
4.1.2 使 用 MacCallum 法 分 析 模 型 拟 合 检验 的 先 验 检验 力 

HALL Satorra-Saris 法 , 使 用 MacCallum 法 分 析 模 型 拟 合 检验 先 验 检验 力 的 过 程 要 简单 
许多 , 只 需要 确定 零 假设 和 备 择 假设 成 立 情况 下 对 应 的 拟 合 指数 、 模 型 自由 度 df、 显著 性 水 
F a 和 统计 检验 力 1-6 五 个 指标 即 可 完成 分 析 。 在 设 定 备 择 假设 对 应 的 拟 合 指数 时 , 研究 者 
可 参考 既往 研究 中 提出 的 关于 该 指数 “大 小 ”的 标准 。 例 如 对 RMSEA 而 言 ， 有 研究 者 指出 ， 
RMSEA 小 于 0.05 为 良好 拟 合 , 0.05 至 0.08 为 合理 拟 合 , 0.08 至 0.1 为 可 接受 的 拟 合 ,大 于 


=a! 
if 


或 


0.1 为 拟 合 较 差 (MacCallum et al., 1996; Yuan & Chan, 2016; Jak et al., 2021)。 对 于 本 文中 给 出 


的 例子 ， 控制 I[ 类 错误 水 平 为 0.05, 若 研 究 者 若 希望 在 RMSEA 的 真 值 为 0.08 时 ， 正 确 拒 绝 


该 模型 的 几率 不 小 于 95%， 此 时 使 用 semPower 包 提 供 的 semPower.aPriori 函数 可 以 计算 得 


出 需要 的 样本 量 为 134 人 。 


4.1.3 使 用 Monte Carlo 法 分 析 模 型 拟 合 检 验 的 先 验 检验 力 


使 用 Monte Carlo 法 分 析 模 型 拟 合 检验 先 验 检验 力 时 确 


同 Satorra-Saris 法 一 致 ,这 里 直接 使 用 4.1.1.1 中 生成 的 到 


定理 论 协 方差 矩阵 的 步骤 和 原则 
E 论 协 方差 矩阵 。 本 例 中 ， 以 等 效 性 


检验 为 分 析 对 象 ， 控 制 I 类 错误 水 平 a 为 0.05, 重复 次 数 R 为 1000。 分 别 取 NW 为 101,201,301 


的 情形 ,计算 在 给 定 N 的 情况 下 1000 次 重复 试验 中 模型 拟 合 指数 落 在 以 RMSEA = 0.08 为 


经 验 临界 值 的 等 效 性 检验 拒绝 域 的 比率 。 需 要 注意 的 是 , 在 Monte Carlo 模拟 的 过 程 中 ,有 


UD 


些 样本 会 导致 结构 方程 模型 在 计算 过 程 中 出 现 迭 代 不 收敛 、 残 差 方 差 为 负数 、 潜 变量 协 方差 
和 矩阵 非 正 定 等 “不 合理 ”的 情况 , 即 出 现 了 不 恰当 解 ,这 些 解 应 在 在 后 续 的 计算 过 程 中 了 予以 


排除 。 本 例 中 , N 分别 为 101 和 201 H 


对 对 应 的 恰当 解 比 例 分 别 为 99.6% 和 100%， 对 应 的 统计 


检验 力 分 别 为 0.724 和 0.996; 对 应 的 各 参数 估计 情况 、 相 对 误差 以 及 95% 置 信 区 间 履 盖 率 
参见 补充 材料 。 本 例 中 全 部 参数 均值 和 标准 误 的 估计 误差 均 不 超过 10%, 95% 置 信 区 间 履 盖 
率 均 处 于 0.91 至 0.98 之 间 , 提示 Monte Carlo 模拟 结果 的 可 信 度 较 高 。 进 一 步 根据 线性 插值 
法 可 求 得 统计 检验 力 为 0.95 时 对 应 的 样本 量 约 为 185 人 。 针 对 同一 个 假设 检验 进行 统计 检 
验 力 分 析 时 , Monte Carlo 方法 得 到 的 结果 看 似 与 Satorra-Saris 法 有 较 大 差异 , 但 这 主要 是 由 


于 线性 插值 带 来 的 误差 。 


事实 上 , 若 通过 缩减 采样 点 间 昌 


E 的 方式 提高 线性 所 


新 求 得 待 估 样 本 量 约 为 151 Ao 


4.2 模型 参数 检验 的 先 验 检验 力 分 析 
4.2.1 使 用 Satorra-Saris 法 分 析 模 型 参数 检验 的 先 验 检验 力 
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参照 Schoemann 等 人 (2017) 的 方式 直接 设 定 理论 协 方差 矩阵 也 并 无 不 可 。 


由 于 既往 研究 中 缺乏 可 以 月 


然 能 够 在 一 定 程度 上 基 导 
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使 用 Satorra—Saris 法 分 析 模 型 参数 检验 先 验 检验 力 的 第 一 步 仍 然 是 设 定 理论 协 方差 矩 
阵 。 诚 然 经 典 的 Satorra-Saris 法 强调 通过 备 择 模型 生成 理论 协 方差 矩阵 , 但 Satorra 和 


Saris(1985) 给 出 的 关于 该 方法 原理 的 讨论 中 并 没有 限 秆 


的 精度 ， 可 重 


论 协 方差 矩阵 的 生成 模式 。 因 此 ， 


日 来 直接 重 构 理 论 协 方差 矩阵 的 数据 ， 而 通过 元 分 析 技 术 虽 


F 现 有 研究 构建 出 理论 协 方差 矩阵 ,但 这 种 方法 需要 耗费 较 多 的 人 


力 物 力 , 单纯 用 来 进行 统计 检验 力 分 析 未 免 有 些 浪费 。 因 此 ,本文 在 进行 理论 协 方差 矩阵 生 
成 时 ,选择 根据 相关 研究 领域 中 “小 、 中 、 大 ”所 对 应 的 经 验 值 进行 设 定 。 此 外 ,由 于 尺度 
变换 不 影响 对 本 例 中 所 关心 的 系数 ( 显 变 量 间 的 路 径 系 数 ) 显 著 性 的 判定 。 因 此 , 在 指定 理论 


协 方差 矩阵 时 可 进一步 将 其 简化 为 指定 相关 和 矩阵 ,参考 Cohen(1988) 的 观点 ,相关 系数 在 0.1、 
0.3、0.5 时 分 别 对 应 “小 入 “中 六 “大 ”的 效应 量 。 在 本 例 中 , 假定 与 控制 变量 有 关 的 相关 
系数 有 “小 ”的 效应 量 , 关注 的 变量 的 相关 系数 有 “中 ”的 效应 量 ， 以 此 设 定理 论 协 方差 矩 
阵 。 以 对 “情绪 创造 力 一 领悟 社会 支持 一 情绪 调节 自我 效能 感 一 压力 后 成 长 ”这 一 路 径 中 介 
效应 的 Wald 检验 为 分 析 对 象 ， 控制 I 类 错误 水 平 a 为 0.05, 求 得 统计 检验 力 为 0.95 时 对 应 的 
样本 量 约 为 795 人 。 
4. 2. 2 使 用 Monte Carlo 法 分 析 模 型 参数 检验 的 先 验 检验 力 

与 之 前 相似 , 这 里 使 用 4.1.1.4 假设 的 理论 协 方差 矩阵 。 本 例 中 ， 以 对 “情绪 创造 力 一 领 
悟 社会 支持 一 情绪 调节 自我 效能 感 一 压力 后 成 长 ”这 一 路 径 中 介 效 应 的 Wald 检验 为 分 析 对 
R, 控制 I 类 错误 水 平 a 为 0.05, 分别 取 N 为 301,401,501, 601 的 情形 进行 模拟 , 得 到 的 结果 
全 部 为 恰当 解 , 对 应 的 统计 检验 力 为 0.613, 0.884, 0.960, 0.985, 对 应 的 各 参数 估计 情况 、 相 
对 误差 以 及 95% 置 信 区 间 履 盖 率 参见 补充 材料 。 本 例 中 全 部 参数 的 估计 误差 均 不 超过 10%, 
关键 参数 的 估计 误差 均 不 超过 5%, 95% 置 信 区 间 履 盖 率 均 处 于 0.91 至 0.98 之 间 ， 提 示 Monte 
Carlo 模拟 结果 的 可 信 度 较 高 。 进 一 步 根据 线性 插值 法 求 得 统计 检验 力 为 0.95 时 对 应 的 样本 
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5. Wit 
5.1 先 验 检 验 力 分 析 时 的 参数 设 定 问 题 

注意 到 , 研究 者 在 使 用 Monte Carlo 法 或 Satorra-Saris 法 进行 先 验 检验 力 分 析 时 ， 均 需 
要 指定 理论 协 方差 矩阵 。 该 步 又 中 , 不 论 研究 者 是 采用 直接 设 定 协 方差 矩阵 的 方式 还 是 通过 
设 定理 论 模型 进而 生成 协 方差 矩阵 的 方式 ,都 需要 设 定 大 量 参数 ， 这些 参 数 的 设 定 会 对 结果 
产生 一 定 影响 , 且 影 响 规律 较为 复杂 。 例 如 , 有 研究 要 对 一 个 两 维度 ,每 个 维度 包含 4 个 题 
目的 量 表 进 行 验 证 性 因子 分 析 。 假 设 Ho 为 两 维度 间 不 存在 相关 , H 成 立 的 情况 下 两 维度 之 
间 存 在 相关 系数 为 0.3 的 正 相关 ， 进行 先 验 检验 力 分 析 时 ， 当 研究 者 设 定 载荷 为 0.7、0.75、 
0.8 时 , 设 定 显 著 性 水 平 为 0.05, 统计 检验 力 为 0.95, 使 用 Satorra-Saris 法 计算 得 到 的 计划 样 
量 分 别 为 527 人、472 人 和 429 A, 使 用 Monte Carlo 法 计算 得 到 的 计划 样本 量 分 别 为 521 
人 、464 人 和 417 Ao XW, 本文 4.2 部 分 给 出 的 例子 中 , 若 将 所 有 控制 变量 与 其 他 变量 的 相 
关系 数 从 0.2 变 为 0.1, 在 其 他 条 件 不 变 的 情况 下 ,使 用 Satorra-Saris 法 计算 得 到 的 计划 样本 
量 为 332 人 , 使 用 Monte Carlo 法 计算 得 到 的 计划 样本 量 为 233 人 。 而 当 分 析 对 象 是 中 介 效 
应 或 有 调节 的 中 介 效应 这 类 由 参数 函数 (通常 是 乘积 ) 构 建 的 指标 时 ,， 若 研究 者 使 用 理论 模型 


生成 协 方差 矩阵 ， 即 使 是 同样 大 小 的 指标 值 ,不 同 的 参数 也 会 导致 分 析 结 果 的 不 同 。 例 如 ， 


对 简单 中 介 效应 模型 而 言 ， 设 定 自 变量 到 因 变 量 的 路 径 系数 为 0.2， 中 介 效 应 值 为 0.16， 显 著 


性 水 平 为 0.05， 统 计 检 验 力 为 0.90, 中 介 前 半路 径 系数 ( 自 变 量 到 中 介 变 量 的 路 径 系数 ) 和 中 


介 后 半路 径 系 数 (中 介 变 量 到 因 变 量 的 路 径 系 数 ) 分 别 取 0.2 和 0.8, 0.4 和 0.4, 0.8 和 0.2( 所 有 


系数 均 为 标准 化 系数 ) 的 情况 下 , 使 用 Satorra-Saris 法 计算 得 到 的 计划 样本 量 分 别 为 258 人 、 
114 人 和 632 A, 使 用 Monte Carlo 法 计算 得 到 的 计划 样本 量 分 别 为 255 人 、88 人 和 620 A. 


从 上 述 例子 中 不 难看 日 


CC 


先 验 检验 力 分 析 中 参数 设 定 对 结果 的 影响 。 这 意味 着 , 研究 者 在 使 用 


Satorra-Saris 法 或 Monte Carlo 法 进行 先 验 检验 力 分 析 时 ， 应 当 审 慎 确 定 参数 的 大 小 ， 必 要 时 ， 


守 的 样本 量 。 
5. 2 理论 协 方差 矩阵 的 产生 方式 


在 缺乏 既往 研究 或 理论 的 支持 下 , 研究 者 可 以 反复 设 定 参数 的 值 ， 并 从 结果 中 选取 一 个 较 保 


如 前 所 述 , 结构 方程 模型 统计 检验 力 分 析 的 Monte Carlo 法 和 Satorra-Saris 法 在 分 析 时 


E 阵 的 指定 存在 两 种 不 同 模式 : 直接 指定 或 通过 指 


需要 指定 理论 协 方差 矩阵 ， 且 理论 协 方差 和 


定 备 择 模型 生成 。 虽 然 从 分 析 过 程 和 分 析 结 果 上 二 者 可 能 并 无 不 同 之 处 , 但 从 逻辑 上 讲 , 不 
同 的 理论 协 方差 矩阵 生成 方式 意味 着 对 结构 方程 建 模 过 程 中 误差 来 源 的 不 同 理解 。Jak 等 人 
(2021) 指 出 作为 衡量 理论 模型 与 真实 数据 间 差 异 的 似 然 比 统计 量 Tu 实际 上 包含 了 两 部 分 的 
误差 : 随机 抽样 造成 的 抽样 误差 和 理论 模型 与 真实 数据 不 匹配 的 误差 。 事 实 上 , 理论 模型 与 


真实 数据 不 匹配 的 误差 还 可 以 再 分 为 两 种 不 同 的 来 源 : 来 自 模 型 内 部 参数 设 定 错误 以 及 模 


型 本 身 对 总 体 的 简化 。 当 研究 者 通过 指定 备 择 模型 的 方式 生成 理论 协 方差 矩阵 时 ， 研 究 者 实 


际 上 默认 了 模型 内 部 参数 设 定 错误 是 导致 误差 的 主要 来 源 , 而 当 研 究 者 通过 直接 指定 的 方 
式 生成 理论 协 方差 矩阵 时 ,， 则 将 这 两 种 不 同 来 源 的 误差 “打包 ”在 了 一 起 。 事 实 上 , 经 典 的 
Satorra-Saris 法 与 MacCallum 法 之 间 也 具有 类 似 的 关系 ,经典 的 Satorra-Saris 法 更 加 强调 来 
自己 经 结构 化 了 的 模型 内 部 的 误差 , MacCallum 法 则 通过 模型 拟 合 指数 将 模型 内 部 的 误差 和 
模型 对 总 体 简化 的 误差 放 在 一 起 考察 ,这 种 细微 的 差别 意味 着 研究 者 需要 在 面 对 不 同 检验 的 


情况 下 灵活 选择 理论 协 方差 矩阵 的 生成 方式 : 当 有 较 强 的 理论 或 既往 研究 指向 误差 可 能 存 
在 于 模型 内 部 时 , 使 用 指定 备 择 模型 的 生成 方式 ; 而 当 研究 者 对 这 些 信 息 掌 握 程度 较 少 时 ， 


使 用 直接 指定 的 方式 或 干脆 换 用 MacCallum 法 进行 分 析 。 


此 外 , 男 一 个 有 趣 的 议题 是 ,使 月 


LU 


先前 的 研究 信息 作为 指定 理论 协 方差 的 依据 是 否 是 合 
适 的 。 诚然 ， 既往 研究 中 确实 存 有 使 用 相似 主题 的 研究 结果 作为 统计 检验 力 分 析 过 程 中 设 定 


参数 依 止 的 先例 , 但 有 研究 者 指出 ,这 种 指定 方式 可 能 错误 地 在 统计 检验 力 分 析 过 程 中 混入 


了 抽样 误差 而 降低 了 统计 检验 力 的 可 靠 性 ,并 进一步 提出 了 纠正 该 误差 的 方案 (Liu & Wang, 


2019)。 对 该 问题 的 深入 讨论 已 超出 了 本 文 的 范畴 ， 有 兴趣 的 读者 可 以 进一步 阅读 Liu 和 


Wang(2019) 的 研究 。 
5.3 非 中 心 和 好 分 布 近似 带 来 的 误差 

Satorra-Saris 法 与 MacCallum 法 这 两 种 分 析 方 法 从 本 质 上 来 说 都 是 利用 非 中 心 ?分布 
作为 当 忆 成 立时 样本 似 然 比 统计 量 Tu 分 布 的 近似 。Steiger Shapiro 和 Browne(1985) 指 出 ， 
这 种 近似 当 且 仅 当 样本 满足 似 然 函数 关于 分 布 的 基本 假设 、 样 本 量 不 太 小 、 拟 合 函数 值 不 是 
太 大 的 情况 下 成 立 。 诚 然 , Satorra, Saris 和 de Pijper(1991) 的 研究 进一步 指出 ， 当 样本 服从 正 
态 分 布 且 理论 模型 与 实际 数据 的 偏离 程度 不 是 十 分 严重 时 ,即使 对 小 样本 而 言 ， 这 种 利用 非 
中 心 x 分 布 进行 近似 的 精度 也 是 基本 令 人 满意 的 。 但 对 于 分 布 假设 不 满足 时 这 种 近似 的 精 
度 是 否 仍 令 人 满意 ， 目 前 仍 缺乏 相关 研究 的 分 析 。 这 意味 着 研究 者 在 对 结构 方程 模型 的 统计 
检验 力 进行 分 析 时 , 应 当 充 分 了 解 样本 分 布 情况 与 所 选 结构 方程 模型 参数 估计 方法 之 间 的 
适 配 性 , 在 进行 统计 检验 力 分 析 和 实际 进行 数据 处 理 时 使 用 同样 的 估计 方法 , 在 样本 分 布 偏 
离 正 态 分 布 的 情况 下 , 也 可 以 考虑 换 用 适用 范围 更 广 的 Monte Carlo 法 进行 统计 检验 力 分 析 。 


前 主流 软件 提供 的 参数 估计 方法 中 , 极 大 似 然 法 (maximum likelihood, ML) 和 广义 最 小 二 


m 


乘法 (generalized least squares, GLS) 要 求 数据 为 正 态 分 布 , 加 权 最 小 二 乘法 (weighted least 


squares, WLS)、 非 加 权 最 小 二 乘法 (unweighted least squares, ULS) 和 对 角 线 加 权 最 小 二 乘法 


(diagonally weighted least squares, DWLS) 对 数据 分 布 要 求 不 高 , 但 ULS 属于 尺度 依赖 方法 


(scale-dependent method), 使 用 时 最 好 保证 各 个 变量 的 测量 单位 统一 , 而 WLS 和 DWLS 通常 
则 需要 1000 以 上 的 样本 量 ( 吴 明 隆 , 2010)。 

此 外 , Satorra-Saris 法 在 使 用 非 中 心 好 分 布 对 样本 似 然 比 统计 量 Tu 的 分 布 进行 近似 时 ， 
同时 使 用 模型 拟 合 “真实 总 体 ” 得 到 的 统计 量 作为 对 非 中 心 参数 的 估计 。 这 意味 着 使 用 
Satorra-Saris 法 分 析 模 型 参数 检验 的 统计 检验 力 时 , 检验 方法 本 身 也 会 对 结果 产生 影响 。 本 
文 4.2 部 分 模型 参数 检验 的 先 验 检验 力 分 析 中 , Monte Carlo 法 与 Satorra-Saris 法 对 同一 检验 
的 分 析 得 出 了 差异 较 大 的 结果 就 是 这 个 原因 。 在 该 例 中 , 在 设 定 a= 0.05, V= 401 的 情况 下 ， 
使 用 Wald 检验 、LM 检验 和 似 然 比 检验 得 到 的 统计 量 对 非 中 心 参 数 进 行 估计 时 得 到 的 非 中 
心 参数 分 别 为 13.399, 22.658 和 23.265, 对 应 的 统计 检验 力 分 别 为 0.955, 0.997 和 0.998。 使 
用 Monte Carlo 法 对 三 种 方法 进行 统计 检验 力 分 析 , 得 到 的 结果 则 分 别 为 0.996, 0.998 和 
0.998。 这 提示 研究 者 在 使 用 Satorra-Saris 法 分 析 Wald 检验 的 统计 检验 力 时 ,可 能 会 由 于 低 


估 统 计 检验 力 而 高 估 所 需要 的 样本 数量 。 


5. 4 Monte Carlo 法 的 精度 与 计算 效率 
诚然 Monte Carlo 法 由 于 可 以 广泛 适用 于 各 类 检验 场景 ， 且 在 模拟 次 数 足 够 多 的 情况 下 
得 出 的 统计 检验 力 的 估计 值 依 概 率 收敛 于 其 理论 真 值 ， 使 得 该 方法 广泛 受到 研究 者 的 青睐 。 


但 Monte Carlo 法 同样 面临 


着 结果 不 稳定 以 及 计算 效率 低下 的 情形 。 与 其 他 两 种 方法 不 同 ， 


Monte Carlo 法 得 到 的 计算 


值 等 于 统计 检验 力 的 真实 值 ， 
可 能 会 得 出 不 同 的 结果 。 


结果 并 不 是 一 个 确定 值 , 而 是 一 个 随机 变量 , 该 随机 变量 的 期 望 


这 就 导致 了 对 同一 个 假定 样本 量 进行 多 次 统计 检验 力 分 析 时 
这 意味 着 研究 者 在 使 用 Monte Carlo 法 进行 统计 检验 力 分 析 时 求 得 


的 样本 量 理应 为 一 约 数 而 非 确切 值 。 换 言 之 , Monte Carlo 法 无 法 像 Satorra-Saris 法 或 


MacCallum 法 那样 给 出 “ 理 #; 


仑 上 达到 目标 统计 检验 力 所 需 的 最 小 样本 量 ”。 这 也 同时 意味 着 


通过 缩短 采样 间距 的 方式 对 基于 Monte Carlo 法 的 统计 检验 力 分 析 的 精度 的 提升 是 有 限 的 。 

此 外 , Monte Carlo 法 的 计算 效率 一 直 以 来 也 是 比较 困扰 研究 者 的 一 个 问题 。 与 Satorra— 
Saris 法 和 MacCallum 法 误差 不 同 , Monte Carlo 法 的 误差 来 源 主要 是 模拟 样本 带 来 的 偶然 误 
差 。 根据 De Moivre-Laplace 中 心 极限 定理 ， 随 着 重复 次 数 R 的 提升 , Monte Carlo 法 所 得 到 


的 统计 检验 力 估 计 值 的 标准 差 会 降低 , 但 其 降低 速率 与 重复 次 数 R 之 间 并 非 线性 关系 , 而 是 


二 次 关系 。 换 言 之 , 若 研 究 者 希望 将 Monte Carlo 法 得 到 的 估计 值 的 精度 提高 一 个 数量 级 ( 相 


当 于 让 估计 值 的 标准 差 降低 为 原来 的 0.1 倍 ),， 则 需要 将 R 扩大 为 原来 的 100 倍 。 这 大 大 降 


低 了 Monte Carlo 法 的 计算 效率 。 并 且 , Monte Carlo 法 本 身 计算 起 来 就 比较 耗 时 , 在 CPU È 
频 为 2.50 GHz 的 笔记 本 电脑 上 运行 本 文中 关于 中 介 模 型 统计 检验 力 分 析 的 R 程序 ， 当 设 定 
HSA R = 1000, 样本 量 N = 401 时 , 使 用 Monte Carlo 法 分 析 Wald 检验 统计 检验 力 所 需 


的 时 间 大 致 为 60s( 单 线程 ); 


所 需 时 间 则 上 升 到 约 900s( 单 线程 ); 若 对 非 参数 Bootstrap 检验 的 统计 检验 力 进行 分 析 ， 即 使 
采用 并 行 计算 技术 将 任务 分 配 到 多 线程 计算 机 的 不 同 线程 上 处 理 ， 所 需 时 间 通 常 也 会 超过 
一 个 小 时 。 此 时 若 反 复 设 定 样本 量 甚 至 参数 进行 多 轮 分 析 从 时 间 上 来 说 是 不 经 济 的 。 诚 然 有 


学 者 指出 此 时 可 以 换 用 参数 


若 将 检验 方法 更 换 为 LM 检验 或 似 然 比 检验 , 执行 同样 的 分 析 


Bootstrap 法 (或 与 之 等 价 的 MC 法 ) 简 化 计算 过 程 ， 节 省 计算 时 


间 (Schoemann et al., 2017), 但 这 种 更 换 也 将 在 一 定 程度 上 造成 统计 检验 力 估 计时 的 误差 。 


5.5 不 同 拟 合 指数 的 选择 对 


MacCallum 法 的 影响 


本 文 主要 介绍 了 用 MacCallum 法 对 RMSEA 相关 的 零 假设 进行 统计 检验 力 分 析 的 理论 
与 实践 。 事 实 上 , MacCallum 法 可 以 用 于 多 种 拟 合 指 数 的 情形 。MacCallum 和 Hong(1997) 将 


该 方法 进一步 扩展 到 了 针对 


是 关于 RMSEA, CFI, MFI 


GFI 和 AGFI 相关 零 假设 的 情形 。Kim(2005) 则 提出 了 当 零 假设 


和 GFI 这 四 种 参数 其 一 时 估算 样本 量 的 简便 公式 。 值 得 注意 的 


是 ， 上 述 参 数 间 的 换算 关系 还 会 受到 模型 本 身 的 影响 。 因 此 ， 目 前 学 界 主流 建议 的 拟 合 指数 
界 值 在 不 同 的 模型 中 并 不 等 价 ,这 就 导致 了 针对 同一 个 计划 模型 ， 采 用 不 同 拟 合 指数 进行 先 
验 检验 力 分 析 时 ， 可 能 得 出 不 同 的 结果 。 例 如 , 根据 温 忠 鹿 等 人 (2004) 的 建议 , 评估 模型 拟 合 
时 , MFI 的 建议 界 值 为 0.85, RMSEA 的 建议 界 值 为 0.08。 此 时 , 对 于 不 同 自 由 度 的 模型 而 言 , 
依据 这 两 个 指标 的 建议 界 值 计算 出 的 计划 样本 量 也 有 所 不 同 。 以 计划 显著 性 水 平 为 0.05, St 
计 检 验 力 为 0.95 的 情况 为 例 ， 当 模型 自由 度 为 10 时 , 根据 MFI 和 RMSEA 计算 出 的 计划 样 
本 分 别 有 77 人 和 383 A; 当 模 型 自由 度 为 50 时, 根据 MFI 和 RMSEA 计算 出 的 计划 样本 分 
别 有 134 人 和 136 A; 而 当 模 型 自由 度 为 100 时 , 根据 MFI 和 RMSEA 计算 出 的 计划 样本 则 
分 别 有 176 人 和 90 人 。 因 此 , 本 文 建议 研究 者 在 使 用 MacCallum 法 进行 先 验 检验 力 分 析 时 ， 
使 用 多 种 拟 合 指数 构建 假设 进行 检验 ,并 在 所 有 计划 样本 量 中 选取 最 大 值 。 至 于 拟 合 指数 的 
有 具体 选择 ,可 以 参考 温 忠 鹿 等 人 (2004) 的 研究 。 
5.6 评价 模型 拟 合 时 的 统计 检验 力 分 析 与 等 效 性 检验 

如 前 所 述 , 既往 研究 中 对 传统 x 检验 的 一 个 批评 就 是 该 检验 将 想 要 验证 的 假设 作为 零 
假设 , 在 结果 报告 中 又 缺乏 对 其 II 类 错误 发 生 概 率 的 分 析 , 使 其 无 法 真正 起 到 验证 模型 的 作 
用 ( 王 阳 等 , 2020). 而 等 效 性 检验 实质 上 就 是 通过 构建 一 个 “可 接受 的 差异 ”标准 的 方式 重 写 
零 假设 ,从 而 将 想 要 验证 的 假设 重新 放 回 了 备 择 假设 的 位 置 , 在 一 定 程度 上 ,也 可 以 将 等 效 
性 检验 视 作 一 种 “颠倒 了 1 与 的 传统 9? Bde”. 根据 MacCallum 等 人 (1996) 的 研究 ， 可 
以 推 知 对 Hoa: Fyio< fo， 且 真实 参数 为 Fyo= ka 的 情形 ， 存在 一 种 与 之 “对 偶 ” 的 情形 : Ho: 
Furo > ka, HIERT AKSO Furo = ko, 这 两 种 情形 下 I 类 错误 概率 和 HI 类 错误 概率 互相 为 
对 方 的 HI 类 错误 概率 与 I 类 错误 概率 。 这 意味 着 , 对 传统 ?检验 而 言 , 在 进行 先 验 检 验 力 分 析 
时 , 若 研究 者 参考 等 效 性 检验 中 的 Ho BFE Fu 的 真 值 ( 即 Fu 的 真 值 不 小 于 eo), 并 设 定 显 著 
性 水 平 为 a, 统计 检验 力 为 1-B 计算 出 计划 样本 量 NM. IA N 同样 可 以 保证 在 控制 显著 性 水 
平 为 B 的 水 平 下 , 对 实际 拟 合 良 好 的 模型 ( 即 Fico = 0), 采用 零 假设 为 Ao: Furo so 的 等 效 性 
检验 进行 检验 时 ,正确 拒绝 零 假 设 的 概率 为 1-a。 因 此 ,本 文 建议 研究 者 在 对 评估 模型 拟 合 
相关 的 假设 检验 进行 统计 检验 力 分 析 时 , 设 定 不 小 于 0.95 的 统计 检验 力 。 这 样 可 以 保障 在 
模型 拟 合 较 差 的 情况 下 , 研究 者 错误 地 认为 模型 拟 合 良好 的 可 能 性 不 高 于 0.05。 


oH 


士 
6. 总 结 


本 文 主要 介绍 了 结构 方程 模型 假设 检验 统计 检验 力 分 析 的 三 种 代表 方法 : Satorra-Saris 
YE. MacCallum 法 和 Monte Carlo 法 , 这 些 方法 各 有 自己 的 适用 情形 。 相 比 于 目前 已 被 广大 


研究 者 所 熟知 的 平均 数 差 异 检验 的 统计 检验 力 分 析 , 结构 方程 模型 的 统计 检验 力 分 析 对 研 
究 者 的 统计 水 平 有 更 高 的 要 求 ， 并且 检 验 过 程 较为 复杂 。 综合 考虑 易 掌 握 和 易 使 用 两 方面 的 
BER, 本 文 建议 研究 者 在 研究 设计 的 过 程 中 根据 自己 需要 进行 的 检验 类 型 选择 恰当 的 统计 
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检验 力 分 析 方 式 ， 具体 可 以 参考 图 2 进行 选择 。 


先 验 检验 力 分 析 还 是 


后 验 检 验 力 分 析 


计划 使 用 的 假设 检验 是 否 为 基于 卡 方 分 布 
方法 的 模型 拟 合 检验 


基于 模拟 或 


重 抽样 的 方法 

优先 推荐 Jonte Carlo 法 

Monte Carlo 法 也 可 参照 先 验 检验 力 分 析 

的 建议 确定 方法 
简单 参数 (路 径 系数 ) 参照 先 验 检验 力 分 析 
是 否 有 明确 备 择 模型 TERREN (his 的 建议 确定 方法 
简单 参数 F 

JacCallun 法 或 计划 对 该 复杂 参数 进 
直接 指定 理论 协 方差 拒 有 Satorra-Saris 法 计划 | 对 该 复杂 参数 进 
阵 的 Satorz a-Saris 法 行 假设 检验 的 方法 


ERS 
Yala 检验 似 然 比 检验 
Skre Seis 


图 2 结构 方程 模型 统计 检验 力 分 析 方 法 选择 建议 


注 : 未 标注 理论 协 方差 矩阵 生成 方法 的 表示 研究 者 可 以 参照 5.2 部 分 的 建议 自行 设 定 


致谢 : LA AMAL RAR, BAERE ZI, RAR X 
稿 过 程 中 提供 的 诸多 帮助 ， 感 谢 两 位 匿名 审 稿 人 对 本 文 提供 的 宝贵 意见 。 
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Power analysis in structural equation modeling: Principles and 


methods 


ZHAI Hongkunl LI Qiang', WEI Xiaowei! 


(‘Department of Social Psychology, Nankai University, Tianjin 300350, China) 


Abstract: Structural equation modeling (SEM) is an important statistical tool in psychology, 
management, and sociology. However, the lack of analyses and reports of statistical power in 
many studies using SEM has reduced the probative power of the study results. There are three 
main types of statistical power analysis methods for SEM: the Satorra—Saris, MacCallum, and 
Monte Carlo methods. The Satorra—Saris method applies to cases where the alternative model is 
clear, analysis object is simple, and the statistic to be examined is based on the y2 distribution. The 
MacCallum method is applicable to the case of y2-based overall model fit tests with unknown 
alternative models. Further, the Monte Carlo method applies to cases where the analysis object is 
complex, or the test is performed using simulation or resampling methods. In practice, researchers 
should clarify the purposes and methods of the test as well as the availability of alternative models 


as a priority. The research method can then be determined on the basis of the above information. 


Key words: structural equation model, statistical power, fitting test, test of model parameters 


表 1 结构 方程 模型 中 的 假设 检验 


零 假设 成 立时 统 
检验 方法 所 处 的 类 别 体 的 检验 方法 应 用 场景 Ho 统计 量 构成 nee 优点 缺点 
Furo=0 
模型 拟 合 的 妃 检验 。 ”检验 模型 拟 合 - (n-1):-Fur 中 心 化 如 分 布 
= 20 
CRR EP pgn Eg n 
效 性 检验 检验 模型 拟 合 。” 王 与 Z 间 的 差距 (n-1): Fu 非 中 心 尼 分 布 
CHa a KEe 
by 在 Fi a EK = = —]). 一 
Na 比较 的 比较 模型 优 几 me — =0 : : (FMLB 中 心 化 站 分 布 
= a uf Xo AE OA = 之 0B MLA PrI E 
N ACRE 模型 参数 的 似 然 比 sg(9)=0 (n-1) (Fut sx 一 inte am 
À MEF 24) 检验 特定 参数 中 心 化 妇 分 布 强 的 假设 ; 对 统计 
p sT 检验 Lo vie = Lo 无 约束 Fut #2) 2 oe 
布 的 方法 、， RIRE- ER 
es RN Lv 8380)=0 g(®) eas 计算 简便 ; 工具 成 熟  ， oaa. 
(经 典 方 Wald 检验 检验 特定 参数 标准 正 态 分 布 制 ; 面 对 较 为 复杂 
~ a Xo 约束 = Lo 无 约束 SEo(9) 区 
© 法 ) Zz 检验 a : 的 参数 函数 〈 如 中 
LM 检验 检验 特定 参数 = 标准 正 态 分 布 介 效应 ) 表现 不 佳 
要 Lo 约束 = Lo 无 约束 SE, 
基于 闪 分 布 a 不 明 , 但 临界 值 
理论 的 检验 拟 合 指数 的 检验 检验 模型 拟 合 。 、_、 给 定 的 拟 合 指数 ”与 中 心 化 闪 分 布 
记分 布 非 中 ax 有 关 
RII annis 8 >a: 不 明 ,但 临界 值 
rie TO RERAMA ”三 与 20 间 的 差距 ”给 定 的 拟 合 指数 。 与 非 中 心 允 分 布 
效 性 检验 
KF (Bx) 有 关 
= Piao=0 给 定 的 拟 合 指数 
检验 模型 拟 合 ”““” BIA A 放宽 了 对 样本 的 假设 ; 
基于 模拟 技 Monte Carlo 也 = Eo (包括 ?统计 量 ) a ee 
Pree ea oe MC 法 不 需 讨论 目标 统计 量 本 ee 
基于 模拟 术 的 方法 法 检验 特定 参数 s(9)=0 而 E HAA 计算 较为 困难 、 计 
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或 重 抽样 ae Siete ~ — 算 耗 时 较 长 、 部 分 
技术 的 方 Fouro=0 给 定 的 拟 合 指数 方法 要 求 使 用 者 有 
检验 模型 拟 合 OY ， BRAN 灵活 度 高 ,对 统计 量 的 人， 、 
法 基于 重 抽样 本 7 了 = Zn (包括 交 统 计量 ) epee 编程 技术 
Pn Bootstrap 法 ”Bootstrap 法 构造 几乎 没有 限制 ; 对 
技术 的 方法 eee g(0) =0 eee res way 
检验 特定 参数 g(0) 经 验 分 布 样本 分 布 要 求 宽松 


Xo 约束 = Lo 无 约束 


《结构 方程 模型 统计 检验 力 分 析 : 原理 与 方法 》 补 充 材 料 


@ 补 充 材 料 A: 案例 部 分 的 R 代码 


# 载 入 程序 包 


if (!require(lavaan)) {install.packages('lavaan') } 


if (!require (mvtnorm) ) {install.packages ('mvtnorm"') } 


if (!require (semPower) ) {install.packages('semPower') } 


4.1 模型 拟 合 检验 的 先 验 检验 力 分 析 
4.1.1 Satorra-Saris 法 


(1) 设 定 用 于 生成 理论 协 方差 矩阵 的 模型 


287 * XL + 0.90 * X2 4.0.75 * X3'° + 0.55 * X4 
209 * X5 + 0.72 * XO £0.58 * X7 0.7L * X8 ft 0.59 * X9 


genMod = ' 

POS ~~ 0.23 * DES 

POS ~~ 0.40 * ANG 

DES ~~ 0.78 * ANG 

POS =~ 0 

DES =~ 0 

ANG =~ 0.73 * X10 + 0.78 * X11 + 0.57 * X12 

XI ~~= 0.24 * XI 

X2 ~~ 0.19 * X2 

X3 ~~ 0.44 * X3 

X4 ~~ 0.70 * X4 

KB: He O05 RS 

X6 ~~ 0.47 * X6 

XT «~~ 0.54 * X7 

X8 ~~ 0.49 * X8 

XY we O52 FOXY 

X10 ~~ 0.45 * X10 

XIL += 0.39 * XLL 

X12 «~ 0.55 * x12" 
(2) 设 定 用 于 分 析 的 模型 

anaMod = ' 

POS =~ X1 + X2 + X3 + X4 

DES =~ X5 + X6 + X7 + X8 + X9 

ANG =~ X10 + X11 + X12' 
(3) 生成 理论 协 方差 矩阵 


rownames (E1) = paste0('X', 1:12) 
(4) 分 析 模 型 


(5) 提取 关键 参数 
f0 = fit@optimSfx * 2 
(6) 用 于 计算 等 效 性 检验 样本 量 的 函数 


setPowerFindEgTestN = function(f0, 


El = cfa(genMod) @implied$cov[[1]] 


fit = cfa(anaMod, sample.cov = El, 


sample.nobs = 101, likelihood = 'Wishart') 


df, 


alpha = 0.05, beta = 0.8, = 0.08) { 


foo = function (N) {pchisq(qchisq(alpha, df, ncp = e*e*df*N), df, f0*N) - beta} 


n = uniroot (foo, c(10, 1e4))S$root 


return (ceiling (n+1)) 
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# (7) 执行 计算 
setPowerFindEgTestN(f0, df = 51, alpha = 0.05, beta = 0.95) 


#4.1.2 MacCallum 法 
semPower.aPriori(effect = 0.08, effect.measure = 'RMSEA', alpha = 0.05, 
power = 0.95, df = 51) -> RES 


summary (RES) 


#4.1.3 Monte Carlo vk 
#(1) 生成 理论 协 方差 矩阵 的 部 分 与 Satorra-Saris 法 相同 ， 这 里 直接 使 用 了 前 面 的 结果 
# (2) 用 于 计算 统计 检验 力 的 程序 

mcCFA = function(N, E, mod, R = 1OOOL) { 


RES = matrix(0, nrow = R, ncol = 62) 
for(i in 1:R){ 
Sample = rmvnorm(N, sigma = E) 


if ('!is.null(colnames (E) ) ) { 


colnames (Sample) = colnames (E) 
jelse{ 
colnames (Sample) = rownames (E) 
} 
RES[i, ] = tryCatch({ 
tempFit = cfa(mod, Sample, std.lv = TRUE, likelihood = 'wishart') 


ind = sum(unlist(parTable(tempFit)[, 14:15])) 
ind = ind / ind 
c(!is.nan (ind), unlist(parTable(tempFit)[, 14:15]), 
as.numeric(fitmeasures(tempFit, 'chisq'))) 
} 
, warning = function(w){c(0, rep(NaN, 61)) }) 
} 
return (RES) 
} 
EG) 利用 程序 计算 统计 检验 力 
set.seed (1234) 
mcCFA(101, El, anaMod) -> RES1 
mcCFA(201, El, anaMod) -> RES2 
#(4) SiR, FRET AE Aad 
resCFA = function(RES, TV, N) { 
R = nrow (RES) 
P = sum(RES[, 1]) / R 


co 


RES = RES[which(RES[, 1] == 1), -1] 

COV = matrix(0, nrow = nrow(RES), ncol = 30) 

for(i in 1l:nrow(RES)) {for(j in 1:30) {COV[i, j] = (TV[j] <= qnorm(0.975, RES[i, jl], 
RES[i, 3+30])) & (TV[j] >= qnorm(0.025, RES[i, j], RES[i, 3+30]))}} 

ME1 = apply(RES[, 1:30], 2, mean) 

ME2 = apply(RES[, 31:60], 2, mean) 


SE = apply(RES[, 1:30], 2, sd) 

PC = apply(COV, 2, mean) 

PM cbind(ME1, ME2, SE, PC, abs(ME1 - TV) / TV, abs(ME2 - SE) / SE) 
RES = RES[, 61] 
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Power = sum(RES <= qchisq(0.05, 51, 0.08*0.08*51* (N-1))) / length (RES) 

colnames (PM) = c('MC MEAN', 'MC_SE', 'MEAN SE', 'Coverage', 'Mean Bias', 'SE Bias") 
RES = list(Proper = P, CHECK = PM, POWER = Power) 

return (RES) 


#(5) 给 出 参数 的 真实 值 
TV = parTable(cfa(anaMod, sample.cov = El, sample.nobs = 101, 
std.lv = TRUE, likelihood = 'Wishart'))[, 14] 

#(6) 给 出 规范 化 输出 
resCFA(RES1, TV, 101) 
resCFA(RES2, TV, 201) 
# (7) 补充 中 点 采样 
mcCFA (151, El, anaMod) -> RES3 
resCFA(RES3, TV, 151) 


4.2 模型 参数 检验 的 先 验 检验 力 分 析 

4.2.1 Satorra-Saris 法 ; Wald 检验 

(1) 设 定理 论 协 方差 矩阵 
E2 = matrix(c( 


EQ Onde. Oly Oecd OL 0.15, 0.01, 
0sty LaO Ok, O22, Only Only aly 
0sty Del, Tad Oly Osdly Orly O21, 
OL, OL 0.0, 120, 023, 0.3;.:0.3; 
Oly Oily Qld Qa 13:0) 063; 043; 
Ost Osby Oh, 0.3) Ves, Lr. Ou3) 
Ody Daly Oe ly Od Oud On3y dad 
), nrow = 7) 

colnames(E2) = c('Age', 'Ser', 'STR', 'ECI', 'PSS', 'RES', 'SRG') 
# (2) 设 定 分 析 模 型 

anaMod = ' 


SRG ~ b*RES + d*PSS + ECI + STR + Ser + Age 

RES ~ k*PSS + ECI + STR + Ser + Age 

PSS ~ a*ECI + e * STR + Ser + Age 

akb := a*k*b!' 

EG) 用 于 计算 统计 检验 力 的 程序 

SSSEM = function(E, mod, PID, alpha = 0.05, beta = 0.95) { 


foo = function (N) { 


fit = sem(mod, sample.cov = E, sample.nobs = N, likelihood = 'Wishart') 
NCP = (parTable (fit) [PID, 14] / parTable (fit) [PID, 15]) ** 2 
return (pchisq(qchisq(1l - alpha, 1), 1, NCP, FALSE) - beta) 


} 
n = uniroot(foo, c(10, 1e4))$root 
return(ceiling(n + 1)) 


} 


# (4) 执行 计算 
ssSEM(E2, anaMod, 29) 
# a i ed i i E eee eae See ee ed i eee Se eS ee SS eS SSS Se wh eS Se eS Se ee Se eS eS Se SS ee Se ed i i 


#4.2.2 Monte Carlo 法 ; Wald 检验 
+O) 生成 理论 协 方差 矩阵 的 部 分 与 Satorra-Saris 法 相同 ， 这 里 直接 使 用 了 前 面 的 结果 
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# (2) 用 于 计算 统计 检验 力 的 程序 


mcSEM = function(N, E, mod, R = 1000L) { 
HVE: 基于 无 完全 共 线 性 的 真实 数据 协 方差 矩阵 的 显 变量 回归 理论 上 不 可 能 出 现 不 恰 
RES = matrix(0, nrow = R, ncol = 58) 
for(i in 1:R){ 
Sample = rmvnorm(N, sigma = E) 
if ('!is.null(colnames (E) )) { 
colnames (Sample) = colnames (E) 


}else{ 


colnames (Sample) = rownames (E) 


} 


解 ， 因 此 没有 考虑 对 其 进行 统计 


for(j in l:ncol(Sample)) {Sample[, j] = (Sample[, j] - mean (Sample[, 
sd(Sample[, j])} 
tempFit <- sem(mod, Sample, likelihood = 'wishart') 
RES[i, ] = unlist(parTable(tempFit)[, 14:15]) 
} 
return (RES) 
} 
EOG) 利用 程序 计算 统计 检验 力 
set.seed (4321) 
mcSEM(601, E2, anaMod) -> RES1 
mcSEM(501, E2, anaMod) -> RES2 
mcSEM(401, E2, anaMod) -> RES3 
mcSEM(301, E2, anaMod) -> RES4 
# (4) 统计 结果 ， 并 进行 规范 化 输出 
resSEM = function(RES, TV) { 
R = nrow(RES) 
COV = matrix(0, nrow = nrow (RES), ncol = 29) 
for(i in 1l:nrow(RES)) {for(j in 1:29){COV[i, j] = (TV[j] <= qnorm(0.975, 
RES[i, j+29])) & (TV[j] >= qnorm(0.025, RES[i, j], RES[i, j3+29]))}} 
El = apply(RES[, 1:29], 2, mean) 
E2 = apply(RES[, 30:58], 2, mean) 
SE = apply(RES[, 1:29], 2, sd) 
PC = apply(COV, 2, mean) 
PM = cbhind(ME1, ME2, SE, PC, abs(ME1 - TV) / TV, abs(ME2 - SE) / SE) 
Power = sum(abs(RES[, 29] / RES[, 58]) >= qnorm(0.975)) / R 
colnames (PM) = c('MC MEAN', 'MC_SE', 'MEAN SE', 'Coverage', 'Mean Bias', 
RES = list (CHECK = PM[1:15, ], POWER = Power) 
return (RES) 
} 
#(5) 给 出 参数 的 真实 值 
TV = parTable(cfa(anaMod, sample.cov = E2, sample.nobs = 101, 
likelihood = 'Wishart'))[, 14] 
# (6) 给 出 规范 化 输出 


resSEM(RES1, TV 


( 
resSEM(RES2, TV 
( 
( 


resSEM(RES3, TV 
resSEM(RES4, TV 
+ 结果 显示 N = 301/601 时 精度 不 足 


) 
) 
) 
) 


RES [i, 


j])) / 


jl, 


"SE Bias") 
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Lim 


#(7) 


mcSEM(601, E2, anaMod) 
mcSEM(301, E2, anaMod) 


resSEM(RES1, 1 


TV) 


resSEM(RES4, 1 


TV) 


新 抽取 精度 不 足 的 样本 


-> RESI 
-> RES4 


@ 补 充 材 料 B: 统计 检验 力 的 后 验 分 析 实 例 


本 部 分 同 


本 部 分 的 演示 样 例 同样 来 自 Zhai 等 人 (2021) 的 研究 。 该 研究 样本 量 为 439 A, 其 中 男 157 A, 女 282 A. 
以 该 研究 中 对 情绪 调节 自我 效能 感 量 表 的 验证 性 因子 分 析 与 以 压力 后 成 长 为 因 变 量 的 链 式 中 介 


模型 为 分 析 对 象 ， 对 其 进行 后 验 统计 检验 力 分 析 。 该 研究 中 ,关于 情绪 调节 自我 效能 感 量 表 的 验证 性 因子 分 
析 结 果 为 : 92 = 151.979, df= 51, RMSEA = 0.067, 参考 王 阳 等 人 (2020) 的 建议 ， 以 RMSEA = 0.08 为 经 验 临界 


值 进行 等 效 性 检验 得 p=0.046。 链 式 中 介 模 型 所 涉及 变量 间 的 协 方差 算 阵 参见 表 B-1， 有 基体 的 模型 设置 同 正 


` 


Xë 


Bat, 路径“ 


绪 创 造 力 一 领 


A 


悟 社会 支持 一 情绪 调节 自我 效能 感 一 压力 后 成 长 ”中 介 效 应 值 为 0.009( 标 准 


化 结果 为 0.021), Wald 检验 z = 2.713, p = 0.007。 参 考 正文 图 2 中 的 建议 , 综合 确定 使 用 MacCallum 法 对 模 


型 拟 合 检验 部 分 i 


行 后 验 检验 力 分 析 , 使 用 Monte Carlo 法 对 链 式 中 介 效 应 部 分 进行 后 验 检验 力 分 析 。 


表 B-1. 链 式 中 介 模 型 所 涉及 变量 间 的 协 方差 矩阵 


1 2 3 4 5 6 7 
1. 年 龄 36.818 
2. 疫情 严重 程度 -0.194 0.936 
3. 疫情 相关 压力 -0.643 2.788 215.412 
4. 情绪 创造 力 -9.107 0.098 9.321 129.108 
5. 领悟 社会 支持 -3.390 -0.992 0.305 60.613 100.388 
6. 情绪 调节 自我 效能 感 -1.217 -0.294 2.198 42.998 29.468 31.922 
7. 压力 后 成 长 -1.981 -0.275 -6.180 24.855 21.677 11.936 21.294 


RMSEA 真 值 为 0.067 时 , 样本 的 ?统计 量 落 在 等 效 性 检验 零 假设 Ho: RMSEA! > 0.08 拒绝 域 中 的 概率 。 根 


1. 模型 拟 合 检验 的 后 验 检验 力 分 析 
在 后 验 检 验 力 分 析 中 , 通常 直接 以 得 出 的 结果 作为 真 值 。 故 本 例 相 当 于 检验 样本 量 为 439、a = 0.05、 


据 正 文中 的 公式 (1a) 和 (2a), 可 以 求 出 后 验 检验 力 , 对 应 的 R 脚本 如 下 : 
(439-1) * 51 * 0.080**2 
(439-1) * 51 * 0.067**2 


ncp0 


ncpa 


Tmlc = 
(POW = 


qchisq(0.05, 
pchisq(Tmlc, 


51, ncp) 
51, ncpa) ) 


运行 上 述 代码 , 可 得 后 验 检验 力 为 0.546。 


验 为 分 析 对 象 ， 控 人 


2. 


介 效 应 检验 的 后 验 检验 力 分 析 


以 对 “情绪 创造 力 一 领悟 社会 支持 一 情绪 调节 自我 效能 感 一 压力 后 成 长 ”这 一 路 径 中 介 效 应 的 Wald 检 


由 类 错误 水 习 


F a 为 0.05, 设 定 N 为 439 进行 分 析 , 得 到 的 结果 全 部 为 恰当 解 ， 对 应 的 统计 


检验 力 为 0.847,， 对 应 的 各 参数 估计 情况 、 相 对 误差 以 及 95% 置 信 区 间 履 盖 率 参见 表 B-2。 本 例 中 关键 参数 
的 估计 误差 均 不 超过 5%, 95% BARKER IAF 0.91 至 0.98 之 间 ， 部 分 参数 的 均值 误差 较 大 ,这 主要 


E 


AE 


jami 


HFSS HH 


RUA) SEN), Sb LIMES BERN RAYNE 5%, 并 且 这 些 参数 都 是 与 控 
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制 变量 相关 的 参数 ,说 明 Monte Carlo 模拟 结果 仍 具 可 信 度 。 这 部 分 的 R 脚本 为 : 
载 入 程序 包 
if (!require(lavaan) ) {install.packages('lavaan') } 
读 入 样本 协 方差 矩阵 
E = matrix (c ( 
36.818, -0.194, -0.643, -9.107, ~-3.390, -1.217, -1.981, 
-0.194, 0.936, 2.788, 0.098, -0.992, -0.294, -0.275, 
-0.643, 2.788, 215.412, 9.321, 0.305, 2.198, -6.180, 
-9.107, 0.098, 9.321, 129.108, 60.613, 42.998, 24.855, 
-3.390, -0.992, 0.305, 60.613, 100.388, 29.468, 21.677, 
-1.217, -0.294, 2.198, 42.998, 29.468, 31.922, 11.936, 
-1.981, -0.275, -6.180, 24.855, 21.677, 11.936, 21.294), 
rownames (E) = c('Age', 'Ser', 'STR', 'ECI', 'PSS', 'RES', 
# 设 定 分 析 模 型 
anaMod = ' 


SRG ~ b*RES + d*PSS + ECI + STR + Ser + Age 


RES ~ k*PSS + ECI + STR + Ser + Age 


PSS ~ a*ECI + e * STR + Ser + Age 
akb := a*k*b! 
TV = parTable(sem(anaMod, sample.cov = E, 


# 执 行 分 析 ， 这 里 利用 wishart 分 布 蔡 代 了 原始 数据 
mod, R = 1000L) { 


sample.nobs = 439, 


mcSEM = function(N, E, 


RES = matrix(0, nrow = R, ncol = 58) 
SaS = rWishart(R, N-1, E) / (N-1) 
for(i in 1:R){ 

tempE = SaS[, , i] 

rownames (tempE) = rownames (E) 


tempFit = sem(mod, sample.nobs = N, 


RES[i, ] = unlist(parTable(tempFit) [, 14:15]) } 
return (RES) } 
set.seed (1234) 
RES = mcSEM(439, E, anaMod) 
resSEM = function(RES, TV) { 
R = nrow (RES) 
COV = matrix(0, nrow = nrow(RES), ncol = 29) 
for(i in l:nrow(RES)) {for(j in 1:29) { 
COV[i, j] = (TV[j] <= qnorm(0.975, RES[i, j], RES[i, 
(TV[j] >= qnorm(0.025, RES[i, j], RES[i, j3+29]))}} 


sample.cov = tempE, 


nrow = 7) 


"SRG') 


likelihood = 


j+29])) 


likelihood = 


& 


'Nishart')y [y 


'wishart') 


14] 


El = apply(RES[, 1:29], 2, mean) 
E2 = apply(RES[, 30:58], 2, mean) 

SE = apply(RES[, 1:29], 2, sd) 

PC = apply(COV, 2, mean) 

PM = cbhind(ME1, ME2, SE, PC, abs(ME1 - TV) / TV, abs(ME2 - SE) / SE) 

Power = sum(abs(RES[, 29] / RES[, 58]) >= qnorm(0.975)) / R 

colnames (PM) = c('MC MEAN', 'MC_SE', 'MEAN SE', 'Coverage', 'Mean Bias', 'SE Bias") 

RES = list (CHECK = PM[1:15, ], POWER = Power) 

return (RES) } 

resSEM(RES, TV) 
表 B-2. 中 介 效 应 Monte Carlo 检验 的 参数 估计 情况 
路 径 理论 值 估计 均值 ”均值 标准 误 。 标准 误 均值 覆盖 率 均值 误差 ” 标准 误 误差 

RES—SRG 0.145 0.144 0.045 0.045 0.949 0.002 0.002 
PSS—SRG 0.118 0.118 0.022 0.022 0.953 0.002 0.023 
ECI—SRG 0.090 0.090 0.023 0.023 0.945 0.003 0.006 
STR 一 SRG -0.034 -0.033 0.013 0.012 0.958 -0.012 0.013 
Ser 一 SRG -0.035 -0.047 0.193 0.195 0.949 -0.324 0.009 
Age 一 SRG -0.017 -0.018 0.030 0.031 0.945 -0.065 0.023 
PSS—RES 0.126 0.127 0.023 0.022 0.961 0.009 0.035 
了 CI 一 人 ES 0.277 0.278 0.020 0.021 0.941 0.002 0.014 
STR 一 人 RES 0.001 0.001 0.013 0.014 0.943 0.574 0.010 
Ser—RES -0.203 -0.201 0.205 0.201 0.950 -0.007 0.018 
Age—RES 0.046 0.046 0.032 0.032 0.949 0.007 0.011 
ECI—PSS 0.472 0.470 0.036 0.036 0.944 0.004 0.008 
STR 一 PSS -0.005 -0.004 0.028 0.029 0.942 -0.126 0.021 
Ser 一 PSS -1.091 -1.091 0.423 0.431 0.941 0.000 0.017 
Age—PSS 0.019 0.023 0.067 0.068 0.946 0.198 0.021 
VE: 加 粗 的 表示 该 路 径 系数 是 关键 参数 ; RES: 情绪 调节 自我 效能 感 , PSS: 领悟 社会 文 持 , ECI: 情绪 创造 力 , SRG: E 


力 后 成 长 , STR: 疫情 相关 压力 , Ser: 疫情 严重 程度 , Age: 年 龄 
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@ 补 充 材 料 C: MacCallum 法 适用 的 拟 合 指数 及 其 与 非 中 心 参数 的 换算 公式 


根据 正文 部 分 的 论述 , 不 难 发 现 MacCallum 法 的 核心 在 于 将 拟 合 指数 根据 其 定义 式 逆 变换 为 非 中 心 参 


数 ， 作 为 对 5 的 估计 值 。 同 时 , 根据 正文 中 的 (3) 式 , MEA: 


d 
ô = (N — 1)Fmo = (N - 1) (Fon — = Ty, — af 


Pa, 如 果 研 究 者 可 以 将 拟 合 指数 “还 原 ” 为 原始 的 ? 值 ， 并 且 知晓 模型 自由 度 ， 那 么 


(C.1) 


理论 上 研究 者 就 可 


以 根据 该 参数 求 出 非 中 心 参数 6 的 估计 值 。 侯 杰 泰 等 人 (2004) 的 书 中 给 出 了 绝 大 多 数 结构 方程 模型 拟 合 指数 
的 定义 ， 从 中 不 难看 出 ， 当 知晓 信息 足够 多 的 情况 下 ,除了 少数 直接 基于 模型 隐 含 协 方差 矩阵 定义 的 拟 合 指 


数 (如 RMR、SRMR 等 ) 外 ,其 余 的 拟 合 指数 均 可 以 “还 原 ” 为 原始 的 入 值 并 进一步 计算 


8 非 中 心 参数 5。 但 


不 同 拟 合 指 数 的 使 用 频率 不 同 , 将 其 换算 为 非 中 心 参数 所 需要 的 额外 参数 (或 者 说 ， 换 算 的 “麻烦 程度 ”) 也 
SE, 故 并 不 是 所 有 的 拟 合 指数 都 在 实践 中 被 应 用 ， 目 前 较为 常用 的 拟 合 指数 有 :RMSEA、MFI、GFI、AGFI 
和 CFI。 这 里 特别 需要 注意 的 是 , 目前 软件 和 论文 中 对 GFI 和 AGFI 的 缩写 存在 一 定 混用 的 情况 ,可 能 不 加 


区 别 地 将 “Goodness of fit index/LISREL GFI”( 下 称 GFID) 和 “Unbiased goodness of fit index/Steiger’s y” (F 


称 GFI“); “Adjusted goodness of fit index”( 下 称 AGFI')#ll “Unbiased adjusted goodness 


of fit index” (FAK 


AGEFID 混 用 。 以 及 为 例 ,R 语言 结构 方程 建 模 程 序 包 lavaan 中 关于 GFI 的 定义 采用 了 GFE, 而 semPower 包 
中 则 采用 了 GFT 。 特 别 需要 说 明 的 是 , MacCallum 法 中 使 用 的 是 GFA AGFI*。 因 此 若 对 lavaan 中 得 到 的 


GFI 值 进行 后 验 检验 力 分 析 ,， 则 可 能 出 现 不 正确 的 结果 。 这 五 种 指数 与 非 中 心 参数 换算 的 公式 为 : 


ô = (N — 1) - df - RMSEA? 
6 = (N — 1) - (—21n MFI) 


E p(1 — GFI*) 
TSWE 2GFI* 


p(1 — AGFI*)df 
p(p +1) — 2df(1 — AGFI*) 


6 = (1 — CFI[(N — 1)Fyuu — dfvuul 


8=(N-1) 


(C. 2a) 


(C. 2b) 
(C. 2c) 


(C. 2d) 


(C. 2e) 


其 中 , p 为 协 方差 矩阵 的 秩 ( 即 变量 维度 ), Fww 和 dfivu 分 别 为 虚 模型 (Null ModeD) 的 拟 合 函数 值 和 自由 度 ， 其 


中 Fyun = 一 In|R|, R 是 变量 的 相关 和 矩阵 。 
参考 文献 : 


ERER, WER, 成 子 娟 . (2004). BYTE RAIA. 北京 : 教育 科学 出 版 社 . (pp. 162-165) 
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@ 补 充 材 料 D: 基于 G*Power 的 结构 方程 模型 统计 检验 力 分 析 

本 文正 文部 分 在 介绍 结构 方程 模型 统计 检验 力 分 析 时 主要 使 用 了 R 语言 作为 软件 工具 。 但 R 语言 本 身 
学 习 成 本 较 高 ,对 于 不 熟悉 编程 的 研究 者 而 言 不 是 非常 友好 。 事 实 上 , R 语言 并 非 进行 结构 方程 模型 统计 检 
验 力 分 析 的 唯一 途径 ,对 于 那些 不 熟悉 R 的 研究 者 ， 只 要 研究 者 本 人 能 较为 熟练 地 使 用 目前 市 面 上 常见 的 


结构 方程 建 模 软 件 (如 LISREL. EQS. AMOS, Mplus 等 )， 再 结合 目前 学 界 常用 的 G*Power, 即 可 完成 基于 
E 分 布 的 结构 方程 模型 统计 检验 力 的 分 析 , 这 在 一 定 程度 上 也 能 够 帮助 研究 者 应 对 一 些 不 是 那么 复杂 的 应 
用 场景 。 

1. 使 用 G*Power 对 结构 方程 模型 统计 检验 力 进行 分 析 的 原理 

正如 正文 中 提 到 的 那样 , 基于 x? 分 布 的 统计 检验 力 分 析 的 核心 包括 两 部 分 : 根据 一 些 方式 估计 出 非 
中 心 参数 ; 巴 根 据 非 中 心 参数 结合 公式 (1a) 或 (1b) 算 出 对 应 情况 下 的 统计 检验 力 。 其 中 的 估计 又 可 分 为 基 
于 协 方差 矩阵 理论 真 值 的 Satorra-Saris 法 和 基于 拟 合 指数 的 MacCallum 法 , 这 部 分 是 需要 研究 者 自己 套用 
公式 计算 或 借助 结构 方程 建 模 软件 进行 计算 的 部 分 ,G*Power 主要 针对 的 是 步 又 @ 中 的 内 容 ， 即 在 已 知 零 假 
设 成 立时 对 应 的 非 中 心 参数 6o、 备 择 假设 成 立时 对 应 的 非 中 心 参数 51、 自 由 度 d、I 类 错误 水 平 a 时 ,根据 
公式 (2a) 或 (2b) 先 计算 出 拒绝 域 的 临界 值 ， 再 根据 公式 (1a) 或 (1b) 计 算出 统计 检验 力 的 过 程 。 

2. G*Power 中 的 Generic x? test 

不 难看 出 步骤 凶 的 核心 是 计算 e 分 布 的 分 位 数 和 单 尾 累积 概率 ,在 G*Power 中 ,相关 功能 可 以 在 
Generic x?test 中 找到 。 打 开 G*Power 后， 先 选 择 检 验 族 (TestFamily) 为 x2 tests， 而 后 便 可 在 检验 方式 (Statistical 


test) 一 栏 中 找到 Generic x? test. G*Power 的 Generic x? test 部 分 主要 提供 了 四 个 工具 :Compromise, Criterion, 
Post hoc 和 Sensitivity。 其 中 后 三 个 工具 与 结构 方程 模型 统计 检验 力 分 析 有 关 , 它们 的 作用 分 别 是 : Criterion: 
给 出 统计 检验 力 1-B、 非 中 心 参 数 (E G*Power 中 用 4 AAMA HE df 求 出 对 应 的 I 类 错误 水 平 a; Post 
hoc: 给 出 非 中 心 参 数 5、 I 类 错误 水 平 a 和 自由 度 df 求 出 对 应 的 统计 检验 力 ]-B; Sensitivity: 给 出 I 类 错误 水 
平 a、 统 计 检验 力 1-B 和 自由 度 df 求 出 对 应 的 非 中 心 参数 。 但 特别 需要 注意 的 是 , G*Power 提供 的 这 些 工 
上 共 都 默认 零 假设 成 立时 对 应 的 非 中 心 参数 60=0， 即 零 假设 成 立时 统计 量 的 分 布 为 中 心 化 闻 分 布 。 在 结构 方 
程 模型 假设 检验 的 经 典 方法 中 ,除去 等 效 性 检验 外 的 其 他 方法 均 能 满足 G*Power 的 默认 条 件 , 因此 可 以 直 
接 借用 G*Power 进行 分 析 。 等 效 性 检验 的 情况 较为 复杂 , 若 等 效 性 检验 的 81 = 0 时 ,可 以 借用 正文 5.6 部 分 
的 结论 , 构建 与 其 对 偶 情 形 , 并 且 将 a 与 B 颠倒 代入 求解 ; 若 等 效 性 检验 的 81 > 0 时 , 则 需 借用 中 心 化 驴 分 
布 为 “ 桥 ”， 建 立 起 两 个 非 中 心 闻 分 布 的 联系 。 

3. 具体 操作 方法 与 示例 

如 前 所 述 , 基于 92 分 布 的 统计 检验 力 分 析 的 第 一 部 分 主要 需要 研究 者 自己 套用 公式 计算 或 借助 结构 方 
程 建 模 软件 进行 计算 , 本文 对 此 步骤 不 进行 详 述 , 仅 简单 介绍 本 步 中 需要 注意 的 几 个 要 点 。@ 本 步 中 的 难点 


二 


在 于 选用 Satorra—Saris 法 时 协 方差 矩阵 理论 真 值 的 生成 , 特别 是 基于 模型 的 理论 协 方差 矩阵 真 值 的 生成 ,这 
需要 研究 者 对 结构 方程 模型 的 基本 原理 有 较 好 的 理解 水 平 。 目 前 市 面 上 大 部 分 结构 方程 建 模 的 教材 都 对 这 
部 分 进行 了 相应 的 介绍 ,系统 学 习 过 结构 方程 建 模 的 研究 者 通常 可 以 在 教材 指导 下 自己 完成 理论 协 方差 矩 
阵 的 生成 工作 .另外 , 生成 理论 协 方差 矩阵 时 需要 计算 矩阵 乘法 ， 这 部 分 可 通过 Excel 的 MMULT 函数 完成 。 

@ 不 同 软 件 计算 多 统计 量 的 方法 略 有 不 同 ,， 对 Mplus m E, x =NMFu, 而 对 LISREL、EQS 和 AMOS 而 言 ,x 
= (N-I)Fu, 这 导致 了 在 使 用 不 同 软件 进行 分 析 时 , 研究 者 可 能 需要 针对 软件 的 计算 方法 对 公式 进行 微调 。 

本 部 分 主要 展示 的 是 第 二 部 分 的 内 容 ， 即 已 知 So af Allo, 61 和 1-B 已 知 其 一 时 如 何 使 用 G*Power 求 出 另 一 
个 的 对 应 值 。 参 考 前 述 对 G*Power 中 的 Generic x? test 的 介绍 , 这 里 主要 分 三 种 情况 给 出 具体 的 操作 方法 与 


示例 : 


(1) io=0,5 >0 

OR 1P, 未知 51( 对 应 先 验 检验 力 分 析 ): 选择 检验 族 (Test Family) x? tests， 检 验方 式 (Statistical test) 
为 Generic x2 test, 检验 力 分 析 类 型 (Type of power analysis) 为 Sensitivity, 设 定 df. a 和 1-B, 求解 出 51, 而 后 
根据 Satorra-Saris 法 得 到 的 Fw (EER MacCallum 法 对 应 的 拟 合 指数 结合 补充 材料 C 部 分 的 公式 推出 结果 。 

例 1: Ue df= 51, 控制 a 水平 为 0.05, 求 使 得 当 RMSEA 的 真实 值 不 小 于 0.08 时 正确 拒绝 零 假设 的 概 
率 不 低 于 95% 的 样本 量 Npian。 

解 : 打开 G*Power, 选择 检验 族 为 x? tests, 检验 方式 为 Generic x? test, 检验 力 分 析 类 型 为 Sensitivity， 设 
定 df=51、a=0.05、1-B=0.95, 求 得 非 中 心 参数 为 43.400, 根据 公式 (C.2a), 求 出 样本 量 为 133.965, 向 上 取 


整 得 Netan = 134 人 。 

OBA 31, AK 1-B( 对 应 后 验 检验 力 分 析 、 使 用 线性 插值 估计 样本 量 时 的 先 验 检验 力 分 析 ): 选择 检验 
WEA x? tests, 检验 方式 为 Generic Xx? test, 检验 力 分 析 类 型 为 Post hoc, 设 定 df. a Al di, 求解 出 1-B。 
例 2: 同 例 1, 使 用 设 定 多 个 N 的 方式 结合 线性 插值 法 估计 Npran。 
解 : 打开 G*Power, 选择 检验 族 为 x tests, 检验 方式 为 Generic 7 test, 检验 力 分 析 类 型 为 Post hoc, WE 
df= 51、a = 0.05、 根 据 公式 (C.2a) 计 算出 入 为 51, 101, 151 时 对 应 的 非 中 心 参数 分 别 为 16.32, 32.64, 48.96, 
分 别 代入 , 求 得 对 应 的 1-B 为 0.566, 0.837, 0.975。 使 用 线性 插值 法 求 出 Npia 约 为 142 Ao 

(2) 80= K(N-71) > 0, 81 =0 

DEA LB, 未 知 8o( 对 应 先 验 检验 力 分 析 ): 根据 正文 5.6 部 分 的 结论 , 构建 与 其 对 偶 情形 一 一 0 = 0, 51 
= k(N-1) > 0, 而 后 交换 a 与 p, 选择 检验 族 为 x tests, 检验 方式 为 Generic x? test, 检验 力 分 析 类 型 为 
Sensitivity, 设 定 df. a 和 1B, 求解 出 S, 进而 补充 材料 C 部 分 的 公式 求 得 N。 

例 3: 已 知 df= 51, 控制 a 水 平 为 0.05， 进 行 等 效 性 检验 , wE RMSEA = 0.08, RE RMSEA 的 真实 


值 为 0 时 正确 拒绝 零 假 设 的 概率 不 低 于 80% 的 样本 量 Npn 


解 : 构建 对 偶 情 形 : 50 = 0, 根据 公式 (C.2a) 有 61 = 0.3264 ° (N-11), EH œ = 0.2, 1-8 = 0.95。 打 开 


G*Power, 选择 检验 族 为 x? tests, 检验 方式 为 Generic x? test, 检验 力 分 析 类 型 为 Sensitivity, € df= S1、u 


=0.2、1-B=0.95, 求 得 非 中 心 参数 为 31.564, 进一步 求 得 Npn 约 为 98 Ao 

QAR So, AK 1-B( 对 应 后 验 检验 力 分 析 、 使 用 线性 插值 估计 样本 量 时 的 先 验 检验 力 分 析 ): 根据 正文 
5.6 部 分 的 结论 , 构建 与 其 对 偶 情形 一 一 80 = 0, 51 = KLVN-7) > 0, 而 后 交换 a 5 p, 选择 检验 族 为 x? tests, 检 
验方 式 为 Generic x test, 检验 力 分 析 类 型 为 Criterion, 设 定 df. a 和 51, 求解 出 1-B。 

例 4: 同 例 2, 使 用 设 定 多 个 N 的 方式 结合 线性 插值 法 估计 Npian。 


解 : 构建 对 偶 情 形 : 50 =0, 根据 公式 (C.2a) 有 61 = 0.3264。(V-7)， 此 时 有 1-B’=0.95. WE 入 分 别 为 51， 


76, 101, 求 得 对 应 的 61 为 16.32, 24.48 和 32.64。 打 开 G*Power 选择 检验 族 为 x? tests, 检验 方式 为 Generic 
L test, 检验 力 分 析 类 型 为 Criterion, 设 定 of = 51, 1- = 0.95, 将 51 的 三 个 值 分 别 代 入 , 求 得 对 应 的 o "为 
0.614, 0.370, 0.180。 使 用 线性 插值 法 求 出 Npi 约 为 99 Ao 


(3) 80 = K(N-1) > 0, 61 = q(N-1) > 0, 60> 61 

此 种 情形 下 问题 过 于 复杂 , 在 此 仅 给 出 一 种 情况 下 的 操作 : 

已 知 51, 未 知 1-B( 对 应 后 验 检验 力 分 析 、 使 用 线性 插值 估计 样本 量 时 的 先 验 检验 力 分 析 ): 注意 到 , 给 定 
自由 度 时 中 心 化 闻 分 布 的 右 尾 概 率 等 价 于 其 临界 值 ， 因此 可 以 以 此 为 “ 桥 ” 连 接 起 两 个 非 中 心 妇 分 布 。 首 
先 选择 检验 族 为 x? tests, 检验 方式 为 Generic x test, 检验 力 分 析 类 型 为 Criterion, 设 定 其 统计 检验 力 为 等 交 
性 检验 的 1-a, df, 8o， 求 出 中 心 化 和 分 布 的 右 尾 概率 , 记 为 y。 然 后 选择 检验 族 为 x? tests, 检验 方式 为 Generic 
X test, 检验 力 分 析 类 型 为 Post hoc, 代入 a = y, 非 中 心 参数 为 61 和 自由 度 , 求 得 对 应 的 统计 检验 力 (此 时 的 
统计 检验 力 为 B 而 非 1-B)。 特 别 需要 注意 的 是 , 这 种 以 中 心 化 多 分 布 为 “ 桥 ” 的 方法 受 限 于 G*Power 对 a 
的 计算 限制 ,G*Power 允许 的 a 范围 是 le-6 到 1-1e-6, 这 也 意味 着 此 种 情形 下 G*Power 的 作用 是 有 限 的 ， 对 
于 此 种 情形 , 若非 中 心 9? 分 布 的 自由 度 足 够 大 时 ,也 可 借用 正 态 分 布 对 非 中 心 y 分 布 进行 近似 从 而 求 出 统 
计 检 验 力 的 近似 值 。 

例 5: 已 知 df= 51, 样本 量 为 301， 设 定 RMSEA = 0.08 进行 等 效 性 检验 , 结果 显示 RMSEA] = 0.043, p 
< 0.05, 求 统计 检验 力 。 

解 : 根据 公式 (C.2a) 计 算出 80 = 97.92, 561 = 28.290。 打 开 G*Power, 选择 检验 族 为 x? tests, 检验 方式 为 


Generic x? test, 检验 力 分 析 类 型 为 Criterion, 设 定 df=51、1-B= 1-0.05 = 0.95, 4 = 97.92, Ka 为 1.014e- 


6。 然 后 选择 检验 力 分 析 类 型 为 Post hoc, 设 定 d=51、4 = 28.290、au = 1.014e-6, 求 得 统计 检验 力 为 0.985。 


例 6: 己 知 df=51, 样本 量 为 439, 设 定 RMSEA = 0.08 进行 等 效 性 检验 , 结果 显示 RMSEA] = 0.067, p 
<0.05, 求 统 计 检 验 力 。 
解 : 根据 公式 (C.2a) 计 算出 So = 142.9632, 61 = 100.2753。 打 开 G*Power, 选择 检验 族 为 Xx tests, 检验 方 
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式 为 Generic Xx? test, 检验 力 分析 类 型 为 Criterion, 设 定 df= 51. 1-B = 1-0.05 = 0.95, A = 142.9632, RIF a 为 
4.064e-012。 注意 到 解 得 的 a 远 超出 了 G*Power 允许 的 范围 , 但 同时 注意 到 本 例 中 df 较 大 ,可 以 考虑 使 用 正 
态 分 布 对 非 中 心 闻 分 布 进行 近似 ， 此 时 读 出 a 为 4.064e-012 时 对 应 的 中 心 化 字 分 布 临界 值 为 152.944， 同时 
计算 非 中心 参 数 为 51、 自由 度 为 51 的 非 中 心 好 分 布 的 均值 和 标准 差 ,根据 非 中 心 闻 分 布 的 相关 结论 有 : 

u = 8+ df = 100.2753 + 51 = 151.275; o = V2df +48 = V2 x 51 + 4 x 100.2753 = 22.430  (D.1) 
以 正 态 分 布 对 非 中 心 好 分 布 进行 近似 ， 可 求 得 均值 为 151.275, 标准 差 为 22.430 的 正 态 分 布 制 点 为 152.944 
时 的 左 尾 概率 约 为 0.530, 与 之 前 B 部 分 得 到 的 理论 值 (0.546) 还 是 较为 相近 的 。 


@ 补 充 材料 E: 案例 中 Monte Carlo 模拟 部 分 各 参数 估计 情况 、 相 对 误差 以 及 95% 置 信 区 间 覆 盖 率 


KE-1 模型 拟 合 检验 先 验 检 验 力 的 Monte Carlo 分 析 


模拟 样本 量 n = 101 模拟 样本 量 n = 201 
路 径 理论 值 ”估计 均值 ”标准 误 均值 ” 标准 误 、 ”估计 均值 ” 标准 误 均值 ” 标准 误 
均值 标准 误 。 ”均值 a 误差 误差 均值 ” 标准 i 均值 a 误差 误差 
POS 一 X1 0.870 0.867 0.083 0.084 0.953 0.004 0.018 0.870 0.059 0.058 0.952 0.000 0.009 
POS 一 X2 0.900 0.896 0.082 0.082 0.946 0.005 0.011 0.901 0.058 0.057 0.941 0.001 0.009 
POS 一 X3 0.750 0.745 0.088 0.089 0.947 0.007 0.008 0.752 0.063 0.062 0.952 0.003 0.013 
POS 一 X4 0.550 0.545 0.096 0.096 0.944 0.009 0.006 0.554 0.068 0.068 0.953 0.007 0.002 
DES 一 X5 0.590 0.589 0.099 0.095 0.956 0.001 0.035 0.590 0.070 0.069 0.958 0.001 0.013 
DES 一 X6 0.720 0.719 0.093 0.096 0.939 0.002 0.030 0.725 0.066 0.068 0.942 0.006 0.029 
DES 一 X7 0.580 0.577 0.091 0.096 0.943 0.005 0.049 0.581 0.065 0.066 0.943 0.002 0.017 
DES 一 X8 0.710 0.702 0.094 0.092 0.949 0.011 0.020 0.711 0.067 0.066 0.953 0.001 0.003 
DES 一 X9 0.590 0.587 0.091 0.092 0.947 0.006 0.019 0.593 0.064 0.066 0.946 0.005 0.025 
ANG — X10 0.730 0.731 0.093 0.093 0.946 0.001 0.006 0.733 0.066 0.066 0.947 0.005 0.010 
ANG > X11 0.780 0.779 0.093 0.094 0.946 0.001 0.013 0.780 0.066 0.069 0.934 0.000 0.040 
ANG > X12 0.570 0.574 0.092 0.095 0.937 0.007 0.032 0.572 0.066 0.066 0.948 0.003 0.012 
XI e X1 0.240 0.237 0.057 0.059 0.932 0.013 0.041 0.240 0.040 0.040 0.955 0.001 0.018 
X2 e X2 0.190 0.189 0.055 0.057 0.948 0.006 0.029 0.188 0.039 0.039 0.951 0.011 0.003 
X3 > X3 0.440 0.437 0.072 0.073 0.932 0.007 0.015 0.438 0.051 0.052 0.944 0.005 0.013 
X4 «> X4 0.700 0.699 0.104 0.104 0.924 0.002 0.001 0.696 0.073 0.074 0.940 0.006 0.015 
X5 e X5 0.650 0.642 0.102 0.106 0.934 0.012 0.032 0.646 0.073 0.075 0.939 0.006 0.038 
X6 «> X6 0.470 0.466 0.085 0.088 0.947 0.009 0.031 0.461 0.060 0.062 0.924 0.019 0.023 
X70 X7 0.540 0.532 0.086 0.088 0.929 0.014 0.017 0.536 0.062 0.063 0.941 0.007 0.029 
X8 «> X8 0.490 0.488 0.087 0.088 0.936 0.003 0.013 0.488 0.062 0.061 0.946 0.004 0.014 
X9 «> X9 0.520 0.516 0.085 0.086 0.937 0.007 0.014 0.519 0.060 0.064 0.933 0.002 0.057 
X10 X10 0.450 0.441 0.086 0.091 0.921 0.021 0.055 0.445 0.062 0.063 0.940 0.011 0.022 
X11l X11 0.390 0.384 0.085 0.087 0.950 0.015 0.018 0.391 0.061 0.062 0.949 0.003 0.012 
X12 e X12 0.550 0.539 0.088 0.091 0.929 0.021 0.031 0.546 0.063 0.062 0.939 0.007 0.014 
POS + DES 0.230 0.227 0.111 0.117 0.935 0.014 0.050 0.230 0.079 0.080 0.937 0.002 0.009 
POS > ANG 0.400 0.398 0.103 0.110 0.929 0.004 0.067 0.401 0.073 0.074 0.953 0.002 0.004 
DES ANG 0.780 0.779 0.071 0.072 0.938 0.001 0.016 0.781 0.050 0.050 0.943 0.002 0.000 


表 E-2 模型 参数 检验 先 验 检验 力 的 Monte Carlo 分 析 


模拟 样本 量 n = 301 模拟 样本 量 n = 401 
路 径 EWE 估计 ae 标准 误 gyz = 但 ee 估计 ae 标准 误 max i 7 pe ‘ 误 
均值 标准 误 均值 误差 误差 均值 标准 误 均值 误差 误差 

了 RES 一 SRG 0.180 0.179 0.057 0.057 0.951 0.006 0.001 0.181 0.049 0.048 0.949 0.001 0.026 
PSS—SRG 0.180 0.183 0.057 0.059 0.942 0.017 0.034 0.180 0.049 0.050 0.948 0.002 0.016 
ECI—SRG 0.180 0.180 0.057 0.056 0.948 0.002 0.016 0.179 0.049 0.050 0.942 0.006 0.013 
STR 一 SRG 0.038 0.038 0.053 0.053 0.949 0.010 0.009 0.040 0.046 0.046 0.950 0.036 0.000 
Ser—SRG 0.038 0.040 0.053 0.053 0.955 0.057 0.001 0.037 0.046 0.046 0.959 0.039 0.012 
Age—SRG 0.038 0.039 0.053 0.054 0.949 0.023 0.019 0.039 0.046 0.046 0.953 0.023 0.008 
PSS—RES 0.220 0.220 0.056 0.056 0.955 0.002 0.011 0.221 0.049 0.047 0.959 0.006 0.027 
ECI—RES 0.220 0.222 0.056 0.056 0.949 0.008 0.005 0.219 0.049 0.046 0.960 0.005 0.058 
STR 一 人 RES 0.047 0.047 0.054 0.057 0.940 0.002 0.049 0.045 0.047 0.048 0.947 0.031 0.018 
Ser—RES 0.047 0.046 0.054 0.054 0.949 0.010 0.006 0.046 0.047 0.048 0.947 0.010 0.017 
Age—RES 0.047 0.045 0.054 0.054 0.955 0.046 0.006 0.048 0.047 0.046 0.955 0.031 0.022 
ECI—PSS 0.282 0.282 0.055 0.054 0.961 0.001 0.017 0.283 0.048 0.046 0.958 0.002 0.036 
STR 一 PSS 0.060 0.060 0.055 0.056 0.949 0.010 0.008 0.062 0.048 0.049 0.951 0.038 0.020 
Ser—PSS 0.060 0.060 0.055 0.056 0.949 0.005 0.003 0.060 0.048 0.048 0.943 0.005 0.010 
Age—PSS 0.060 0.055 0.055 0.059 0.935 0.077 0.062 0.058 0.048 0.048 0.942 0.025 0.003 


E: n 为 501 和 601 时 的 参数 估计 情况 见 下 页 续 表 ， 加 粗 的 表示 该 路 径 系数 是 关键 参数 ; RES: 情绪 调 
SRG: 压力 后 成 长 , STR: 疫情 相关 压力 , Ser: 疫情 严重 程度 ,Age: 年 龄 


我 效能 感 , PSS: 领悟 社会 支持 , ECI: 情绪 创造 力 ， 


dt 


表 E-2( 续 ) 模型 参数 检验 9 


E 验 检验 力 的 Monte Carlo 分 析 ( 续 表 ) 


模拟 样本 量 n = 501 模拟 样本 量 n = 601 
路 径 理论 值 估计 均值 ”标准 误 y 均值 标准 误 估计 均值 估计 均值 均值 ” 标准 误 
均值 标准 误 、 ”均值 sara RE RE 均值 ”标准 误 均值 ”标准 误 误差 误差 
RES—SRG 0.180 0.179 0.044 0.043 0.952 0.008 0.013 0.178 0.040 0.041 0.940 0.011 0.022 
PSS 一 SRG 0.180 0.180 0.044 0.046 0.943 0.003 0.035 0.182 0.040 0.041 0.944 0.009 0.023 
ECI-SRG 0.180 0.183 0.044 0.043 0.952 0.015 0.019 0.180 0.040 0.040 0.941 0.001 0.000 
STR 一 SRG 0.038 0.038 0.041 0.040 0.956 0.017 0.019 0.039 0.038 0.038 0.951 0.009 0.006 
Ser—»SRG 0.038 0.039 0.041 0.041 0.954 0.028 0.002 0.038 0.038 0.039 0.947 0.005 0.031 
Age—SRG 0.038 0.041 0.041 0.041 0.957 0.064 0.013 0.039 0.038 0.037 0.947 0.015 0.009 
PSS—RES 0.220 0.221 0.044 0.042 0.962 0.003 0.037 0.218 0.040 0.041 0.951 0.010 0.018 
ECI-RES 0.220 0.221 0.044 0.042 0.958 0.006 0.045 0.220 0.040 0.039 0.952 0.002 0.029 
STR 一 RES 0.047 0.050 0.042 0.043 0.947 0.067 0.020 0.047 0.038 0.038 0.950 0.005 0.019 
Ser—>RES 0.047 0.045 0.042 0.041 0.958 0.029 0.010 0.046 0.038 0.039 0.943 0.017 0.009 
Age 一 RES 0.047 0.047 0.042 0.042 0.947 0.013 0.003 0.048 0.038 0.038 0.955 0.038 0.005 
ECI—PSS 0.282 0.282 0.043 0.042 0.948 0.001 0.010 0.280 0.039 0.037 0.963 0.008 0.049 
STR 一 PSS 0.060 0.059 0.043 0.042 0.956 0.018 0.023 0.061 0.039 0.040 0.947 0.012 0.011 
Ser—PSS 0.060 0.059 0.043 0.043 0.946 0.010 0.010 0.061 0.039 0.041 0.945 0.012 0.041 
Age 一 PSS 0.060 0.060 0.043 0.043 0.945 0.000 0.000 0.059 0.039 0.040 0.948 0.009 0.029 


